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I.  BACKGROUND 


The  testing  and  eharaeterization  of  U.S.  Navy  underwater  acoustie  devices 
remains  a  fundamentally  important  activity  linked  directly  to  both  fleet  and  mission 
readiness.  Testing  of  new  prototypes,  sampling  of  production  line  devices  and  re¬ 
qualification  of  in-service  assets  are  all  components  of  this  readiness.  These  testing 
efforts  represent  measures  put  in  place  to  help  assure  that  the  required  level  of  acoustic 
performance  is  being  achieved. 

The  Underwater  Sound  Reference  Division  (USRD)  represents  the  Navy's 
primary  activity  for  underwater  acoustic  calibration,  test,  and  evaluation  measurements. 
In  this  capacity,  the  USRD  maintains  specialized  measurement  facilities  with  the 
capability  to  simulate  real-world  ocean  environments.  The  USRD  is  part  of  the 
Submarine  Sonar  Department  of  the  Naval  Undersea  Warfare  Center  Division  Newport 
located  in  Newport,  Rhode  Island  (www.npt.nuwc.navy.mil/USRD).^*^ 

Within  the  USRD  is  an  acoustic  measurement  facility  designated  as  the  Acoustic 
Pressure  Tank  Facility  (APTF).  The  APTF  test  vessel  consists  of  a  closed  steel  tank  3.81 
meters  in  diameter  and  11.1  meters  in  length  with  two  access  ports.  The  facility  allows 
the  acoustic  test  and  evaluation  of  devices  and  materials  over  a  wide  range  of  simulated 
ocean  environments.  Temperatures  from  2  degrees  Celsius  to  35  degrees  Celsius  and 
hydrostatic  pressures  up  to  18.6  megapascal  (a  depth  equivalent  of  1860  meters)  can  be 
achieved. 


Naval  Undersea  Warfare  Center  Division  Newport  -  Aconstic  Pressnre  Tank  Facility 
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The  APTF  represents  a  significant  U.S.  Navy  capital  investment.  As  Navy 
demands  and  needs  have  changed  over  the  years,  it  has  become  increasingly  more 
important  to  revitalize  existing  Navy  assets  to  meet  new  requirements.  While  it  is  a 
common  need  to  characterize  the  performance  of  acoustic  devices  under  realistic  ocean 
environmental  conditions,  in  some  instances  these  devices  may  require  test  distances  in 
excess  of  that  achievable  within  the  confines  of  existing  Navy  test  vessels.  As  an 
example,  these  devices  may  entail  designs  that  have  large  apertures  yet  may  also  function 
at  relatively  high  frequencies  (as  the  U.S.  Navy  High  Frequency  Sail  Array).  This 
distance  required  for  proper  acoustic  characterization  is  designated  as  that  distance  where 
far  field  conditions  occur.  Test  distances  that  meet  or  exceed  this  far  field  criterion 
allow  the  device  to  function  as  it  would  in  the  free-field  condition  of  the  open  ocean. 
Testing  at  distances  that  fall  short  of  this  criterion  are  designated  as  being  tested  in  the 
near  field.  Testing  in  the  near  field  is  subject  to  contaminations  that  obscure  the  true  far 
field  performance  of  the  device  and  thus  is  typically  not  desirable. 

This  thesis  effort  addresses  a  technique  that  uses  near  field  acoustic  data  with  a 
suitable  math  model  to  extrapolate  far  field  performance.  The  use  of  “CHIEF  for 
MATLAB”'-^^  coupled  with  a  least  squares  solution  methodology  are  utilized  in  this 
technique.  The  data  output  rendered  from  the  technique  mimics  data  expected  within  test 
vessels  that  would  be  prohibitively  expensive  to  fabricate  and  maintain. 

The  remainder  of  this  thesis  is  organized  as  follows: 

Chapter  II  addresses  the  definition  and  meaning  of  near-field  and  far-  field  conditions. 

Chapter  III  provides  an  overview  of  the  Combined  Helmholtz  Integral  Equation 
Eormulation  (CHIEE)  and  the  A  and  B  matrices  that  are  used  in  subsequent  computations 
in  the  thesis  methodology. 

Chapter  IV  details  the  steps  of  the  thesis  near-field  to  far-field  methodology  and  criteria 
for  assessing  the  data  output. 
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Chapter  V  provides  an  example  solution  set  using  the  methodology. 


Chapter  VI  provides  an  error  analysis  of  the  technique  when  positional  uncertainties  of 
the  various  acoustic  devices  used  during  a  test  exist. 

Chapter  VII  looks  at  other  specific  test  examples  inclusive  of  testing  in  more  extreme 
near-field  cases  and  test  results  expected  from  what  is  likely  to  be  considered  as  a  typical 
test  set-up  and  scenario. 
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II,  NEAR-FIELD/FAR-FIELD  OVERVIEW 


A,  INTRODUCTION 

A  fundamental  physical  dimension  addressed  in  conducting  acoustic 
measurements  is  separation  distance.  The  desired  separation  distance  between  devices  in 
a  test  setup  places  them  at  distances  where  1/r  spreading  occurs  (with  r  defined  as  that 
distance  between  the  projecting  and  receiving  devices).  This  mimics  expected 
performance  in  a  free-field  state  and  enables  an  accurate  characterization  of  the  device 
for  the  typical  ocean  environment  for  which  it  is  designed.  1/r  spreading  gives  the 
expected  6  dB  loss  when  separation  distance  is  doubled  This  type  of  spreading  occurs 
when  “far-field”  conditions  have  been  met. 

Testing  of  a  transducer  with  finite  dimensions  (i.e.,  not  a  point  source)  at 
distances  inside  the  far-field  limit  results  in  performance  characterization  that  is  obscured 
by  phase  superposition  of  the  comprising  elements.  That  is,  the  complex  response 
measured  at  a  particular  point  is  made  up  of  a  number  of  individual  responses  each  with 
widely  varying  phase  terms.  Thus,  the  “near  field”  of  a  particular  set-up  is  characterized 
by  radical  fluctuations  in  pressure  with  changes  in  range  (in  a  given  direction).  This 
obscures  the  nature  of  the  device’s  true  far  field  performance. 

The  following  curve  demonstrates  the  variation  in  the  near  field  and  far  field 
responses  seen  for  a  particular  acoustic  device  (a  linear  array).  The  response  curve  was 
obtained  by  summing  up  the  complex  responses  from  all  source  elements  at  each  of 
several  specified  ranges.  Note  the  widely  changing  response  in  the  near  field  of  the 
device  and  the  1/r  spreading  resulting  at  the  farther  ranges. 
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Figure  1  -  Pressure  Response  with  Change  in  Range 


The  field  data  at  varying  ranges  in  the  previous  figure  (Figure  1)  was  computed 
for  a  30.5  cm  linear  array  using  a  test  frequency  of  600  kHz.  A  far  field  distance  of 
approximately  37  meters  is  estimated  for  this  device.  The  data  set  was  normalized  to  the 
maximum  response  computed.  Computations  were  based  upon  modeling  the  source  as  a 
“Continuous  Line  Source”  as  detailed  in  Kinsler  and  Frey,  Section  7.3  (4'*'  Edition) 

B.  THEORY 

The  properties  of  the  corresponding  acoustic  field  for  a  given  device  can  typically 
be  estimated  from  the  geometry  of  the  source  and  its  wavelength,  X  ,  of  operation.  The 
extent  of  the  near  field  can  be  estimated  through  the  examination  of  the  acoustic 
contributions  from  the  device’s  various  geometric  components.  The  figure  that  follows  is 
used  to  examine  these  details. 
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Figure  2  -  Geometry  used  in  estimating  extent  of  Near  Field 


In  the  above  figure,  ,  represents  the  distanee  from  the  furthest  element  of  the  souree 

(in  this  example  this  is  a  line  array  of  length  L)  to  a  field  point.  The  magnitude,  , 
represents  the  distance  from  the  nearest  element  of  the  source  to  the  same  field  point. 
The  difference  between  these  two  magnitudes  is  defined  as  A  r  =  .  At  great 

distances,  this  difference  (Ar  )  is  minimal.  At  closer  distances  and  as  the  field  point 

approaches  the  source  on  axis,  this  difference  can  be  significant.  When  the  difference 
approaches  XU  (over  the  difference  obtained  at  great  range),  the  phases  of  the  signals 
from  the  individual  points  on  the  source  combining  at  the  field  point  will  be  of 
sufficiently  different  value.  The  net  effect  of  this  is  that  the  axial  pressure.  Pax,  can 
change  radically  over  small  changes  in  distance  (see  Figure  1).'-"^’^*’^ 

Referencing  Figure  2  again  gives  the  following  expression  for  indicating  a 
reasonable  boundary  for  expected  far  field  performance: 

r_/L  =  L/4A  (1) 
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L/2  represents  half  of  the  maximum  dimension  of  the  souree. 


A  rule  of  thumb  for  far  field  separation  distanee  that  has  been  used  in  the  USRD  faeilities 
has  been  based  upon 

D  " 

with  D  equal  to  the  maximum  dimension  of  the  aperture 
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Ill,  CHIEF  -  INTRODUCTION 


A,  BACKGROUND 

The  CHIEF  (Combined  Helmholtz  Integral  Equation  Formulation)  program  was 
originally  developed  at  the  Naval  Ocean  Systems  Center  (NOSC)  from  1965-1970  by 
Drs.  H.A.  Schenck  and  G.W.  Benthien.  It  was  originally  developed  to  compute  acoustic 
radiation  that  would  be  produced  from  a  given  arbitrarily  defined  body.  Over  the  years 
its  capabilities  have  been  expanded  by  others  to  utilize  more  powerful  and  available 
computing  hardware  and  to  include  plane-wave  scattering  and  finite-element  defined 
input  descriptions  . 

The  version  of  CHIEF  used  in  this  effort  was  written  for  implementation  with 
MATEAB  by  Mr.  Stephen  Forsythe,  NUWCDIVNPT,  during  fiscal  years  1999-2000. 
This  represents  a  unique  implementation  of  CHIEF  that  takes  advantage  of  the  matrix 
handling  capabilities  of  a  commercial  off  the  shelf  product.  While  achieving  the 
functionality  of  the  original  CHIEF  program,  Mr.  Forsythe  further  augmented  this 
version  to  utilize  spatially  defined  material  properties  rather  than  just  relying  on  boundary 
value  problem  descriptors  .  The  open  nature  of  the  coding  that  was  implemented  also 
allows  access  to  non-traditional  CHIEF  outputs.  It  is  this  latter  implementation  detail  that 
made  this  version  of  CHIEF  the  proper  choice  for  this  effort.  Specifically,  the  extraction 
of  the  CHIEF  A  and  B  matrices  for  use  outside  of  the  protocols  of  CHIEF  was  the  unique 
feature  that  this  version  provided  and  that  this  thesis  methodology  required. 

B,  CHIEF  THEORY 

This  section  provides  a  brief  glimpse  at  the  definition  and  formulation  of  the 
CHIEF  A  and  B  matrices  and  the  use  of  the  Helmholtz  Integral  formula  to  deduce 
pressure  at  a  point  in  space.  This  thesis  effort  does  not  concentrate  on  the  exacting 
derivation  of  these  parameters  but  rather  provides  an  overview.  More  complete  details 
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may  be  found  on  this  derivation  in  the  doeumentation  provided  with  the  CHIEF  1988 
release  . 


For  both  aeoustie  radiation  and  aeoustie  seattering  problems,  Green’s  theorem 
speeifies  that  the  total  aeoustie  pressure  p  (  pinc  +  Ps,  with  pinc  defined  as  the  ineident 
pressure  on  a  surfaee  S  due  to  sourees  in  the  region  of  interest  and  ps  defined  as  the 
eomplex  seattered  pressure  exterior  to  S)  and  the  normal  veloeity  v  satisfy  the  Helmholtz 
integral  formulas: 


p{x)-p  inc  (x)  ;  X  exterior  to  5* 
^  p{x)  -  pinc{x) ;  X  on 
—  pinc  (x) ;  X  interior  to  S 


with: 


1 

G(x,x’)  =  - - ^ 

4;r  x-x 


(2) 


where: 

G(x,x')  =  three  dimensional  free-spaee  Green’s  funetion  and 

k :  aeoustie  wave  number 
co:  radian  frequeney 
e:  sound  speed  in  the  fluid 

10 


- :  signifies  differentiation  in  the  outward  normal  direetion 

dn 


and  a  harmonie  time  dependenee  of  e'®  is  assumed. 

If  ^  is  some  point  loeated  on  the  surfaee,  equation  (1)  from  above  leads  to: 

,cr)ds{cr)  =  PinX<^)  +  i(op\  v{<j)G{(; ,cr)ds{<j)  (3) 

s  s 

C.  CHIEF  -  NUMERICAL  APPROACH 

Sinee  the  goal  in  CHIEF  is  to  solve  this  problem  numerieally,  the  surfaee  S  is 
broken  up  into  an  Ns  number  of  non-overlapping  “patehes”.  The  pressure  and  veloeity  at 
eaeh  of  these  patehes  is  approximated  by  average  values,  designated  as  Pn  and  Vn, 
respeetively. 

The  diseretization  of  equation  (3)  leads  to: 

1  v,  d 

tr -  S ^ .  or)ds{cj)  =  ) 

2  „=i  on^ 

+  fJ  G(C  ,  cj)ds{cj)  (4) 

«=i  s„ 


The  matriees  A  and  B  are  subsequently  defined  as  follows: 
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A 


mn 


2 


d 

dn^ 


B 


mn 


(5) 


If  the  variables  p  ,  v,  and  pi„c  are  defined  as  vectors  of  surface  pressures, 
velocities  and  incident  surface  pressures,  then  equation  (4)  can  be  more  simply  stated  as: 

A  p  =  B  V  +  Pine  (6) 

Note:  further  refinement  of  the  A  and  B  matrix  definitions  in  equation  (5)  above 
are  required  to  handle  specific  wave  numbers  and  surface  pressure  conditions.  The 
matrix  A  is  otherwise  not  suited  to  achieving  an  over  determined  system  solvable  via 
least  squares  solution  technique.  This  further  refinement  is  not  addressed  here. 

As  mentioned  in  the  early  part  of  this  section,  the  Helmholtz  Integral  relations  can 
also  be  used  to  deduce  pressures  at  defined  locations  in  space.  The  pressure  at  any  field 
point  X  (near  field  or  far  field)  can  be  expressed  in  terms  of  the  surface  pressure  and 
normal  velocity  per  the  following  relationship: 


Psix)  =  Pix)-pinc{x) 


(V) 


G(x,  a)  +  i(i>pv{<j)G{x,  a)  \ds{<j) 


In  like  fashion  as  the  formulation  of  the  CHIEF  A  and  B  matrices,  it  desired  to 
solve  equation  (7)  numerically  by  breaking  the  surface  S  into  a  mosaic  of  non¬ 
overlapping  patches.  This  leads  to  the  expression  identified  in  equation  8  that  follows. 
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Psi^)  =  Pi^)- 


(8) 


"  r  5  "  r 

=  — G(x,cr)J5(cr)  +  z®j!?^F„  jG(x,cr)J5(cr) 

H=1  n=l 

Using  the  condition  where  pinc{x)  is  equal  to  zero  (indicating  that  there  is  no 
incident  pressure  on  the  surface),  the  discrete  notation  of  equation  8  may  be  further 
simplified  to: 

p{x)  =  p  delG  +  V  G  (9) 

where  each  of  the  constituents  (from  left  to  right)  on  the  rhs  of  equation  9  are  vectors  of 
the  pressure  at  a  known  surface  location,  the  integration  of  the  normal  derivative  of  the 
Greens  function  for  the  given  set-up,  the  normal  velocity  at  a  known  surface  location  and 
the  integration  of  the  Greens  function  for  the  given  set-up,  respectively.  Please  note  that 
X  designates  the  three  dimensional  location  in  space  of  a  defined  field  point  and  cr 
(equation  8)  designates  the  three-dimensional  location  in  space  of  the  surface  patches. 

D,  SUMMARY 

CHIEF  utilizes  the  Helmholtz  integral  relations  of  surface  pressures  and  normal 
velocities  along  with  a  defined  surface  incident  pressure  condition.  Computations  within 
CHIEF  proceed  numerically  via  the  re-structuring  of  the  continuous  surface  into  a  mosaic 
of  individual  non-overlapping  surface  patches.  This  re-structuring  allows  the  formulation 
of  the  generalized  equation: 

Ap  =  Bv  +  pinc  (10) 
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with  p  ,  V,  and  pinc  defined  as  veetors  of  surfaee  pressures,  normal  veloeities  and  ineident 
surfaee  pressures,  respeetively.  A  is  the  matrix  formed  from  surfaee  integrals  of  the 
normal  derivative  of  pressure  over  the  surfaee  patehes  and  B  is  the  matrix  formed  from 
surfaee  integrals  of  the  pressure  over  the  surfaee  patehes.  Both  the  A  and  B  matrices  can 
be  formed  from  known  acoustic  constants  and  a  defined  geometry  . 

CHIEF  also  yields  a  mechanism  to  calculate  pressure  at  a  defined  field  point  via 
the  use  of  known  surface  pressures  and  normal  velocities  and  integrated  three- 
dimensional  Greens  function  derived  relationships  (between  surface  locations  and  the 
defined  field  point  x).  This  leads  to  the  following  relationship: 

;7(x)  =  p  delG  +  V  G  (11) 

with  p{x)  defined  as  the  pressure  at  a  designated  field  point  x,  p  and  v  defined  as  vectors 
of  surface  pressures  and  normal  velocities,  respectively,  and  delG  and  G  defined  as  the 
integrated  Greens  functions  relationships  between  field  point  x  and  the  surface  patch 
locations. 

Equations  10  and  11  in  the  above  summary  form  most  of  the  basis  that  will  be 
utilized  in  setting  up  the  near  field-far  field  extrapolation  problem’s  defining  equations. 
The  methodology  will  require  that  Equation  10  be  equated  to  zero  (indicating  that  no 
external  sources  are  present)  and  equation  (11)  will  be  used  to  relate  to  pressures  at 
defined  near  field  points. 

Due  to  the  numerical  approach  that  CHIEF  utilizes,  it  is  important  that  the  number 
of  patches  used  in  the  mosaic  to  cover  the  surface,  S,  be  properly  sized.  A  single 
complex  value  of  pressure  and  velocity  per  patch  location  is  a  required  condition  in 
CHIEF;  these  values  of  pressure  and  velocity  over  the  entire  patch  are  represented  by 
values  computed  for  the  patch  center. 
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IV.  STEPS  OF  THE  METHODOLOGY 


A.  INTRODUCTION 

This  section  outlines  and  details  the  steps  of  the  methodology  required  to  take 
near  field  array  element  data  through  the  processing  steps  required  to  achieve  an 
extrapolated  far  field  response.  The  steps  provide  the  procedural  details  as  well  as  criteria 
for  parameter  selection.  Practical  considerations  of  limitations  of  the  processes  are  also  a 
necessary  part  of  these  steps  since  the  effort  entails  the  implementation  of  the  technique 
on  available  computing  platforms.  This  section  is  supplemented  by  the  “Explanation  of 
MATLAB  Coding”  section  that  appears  in  Appendix  C. 
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CHIEF  Enclosing  Sphere  placed  around  Linear  Transmitting  Array 
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Figure  3  -  Typical  Methodology  Set-up  showing  device  under  test,  receive  array  and 
CHIEF  Enclosing  Sphere 
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The  figure  above  shows  a  typieal  set-up  that  forms  the  basis  of  the  methodology. 
The  transmitting  array  and  the  reeeive  array  both  represent  actual  physical  devices,  while 
the  CHIEF  enclosing  sphere  represents  an  artificial  surface  in  space.  The  methodology 
seeks  to  define  the  pressure  and  velocity  characteristics  of  this  spherical  surface  for 
eventual  extrapolation  to  far  field  pressure. 

While  some  of  the  early  steps  are  not  execution  order  specific,  it  is  practical  to 
treat  all  of  these  steps  as  sequential  in  the  order  that  they  are  presented. 

B,  DEFINE  CHIEF  ENCLOSING  SPHERE 

This  step  of  the  process  takes  into  account  the  physical  size  of  the  acoustic  device, 
and  its  frequencies  of  operation.  Since  the  enclosing  sphere  (formulated  using  “CHIEF 
for  MATEAB”  protocols  )  will  be  formed  from  multiple  repeated  copies  of  a  single 
“slice”  of  a  sphere  (imagine  a  segment  from  a  section  of  a  tangerine),  an  estimate  of  the 
number  of  comprising  slices  required,  ie:  its  symmetry  order  (sph  sym),  is  performed  in 
this  step.  The  number  of  divisions  within  this  single  “slice”  of  the  sphere  (designated  as 
patches)  are  also  formed  in  this  step.  The  slice  will,  therefore,  be  represented  by  a 
sph_numv  x  1  matrix  with  sph_numv  representing  the  number  of  equally  spaced  patches 
in  the  direction  perpendicular  to  the  symmetry  order  .  (Eimitations  are  imposed  on  both 
of  these  quantities  based  upon  computational  capabilities  and  efficiencies  of  both  the 
MATEAB  coding  and  the  standard  COTS  computer  hardware  used  for  this  effort.) 

C.  DEFINE  RECEIVE  ARRAY  AND  PROXIMITY  TO  ACOUSTIC  DEVICE 

This  step  defines  the  number  of  receive  array  elements  (numb)  and  the  array’s 
distance  from  the  acoustic  device  under  test.  It  is  assumed  that  the  elements  in  this  array 
are  in  a  vertical  line  and  are  equally  spaced.  The  technique  does  not  necessarily  require 
this  for  any  particular  technical  reason,  rather  this  is  based  more  upon  the  expected 
configuration  that  will  be  suited  for  an  experimental  set-up.  Maximum  spacing  of  the 
elements  in  this  array  is  expected  to  be  dependent  upon  the  required  frequencies  of 
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operation.  The  maximum  number  of  receive  elements  for  this  technique  will  be 
established  more  by  logistical  considerations  rather  than  for  technical  reasons  (ie:  there 
will  be  some  practical  limit  to  the  number  of  elements  that  can  be  made  available  for  the 
receive  array  based  upon  available  sensors  in  inventory  and  the  number  of  channels 
available  on  the  data  acquisition  system).  The  distance  of  the  receive  array  from  the 
acoustic  device  will  be  in  the  near  field  of  that  device.  Part  of  the  evaluation  of  this 
technique  will  be  to  examine  final  extrapolated  data  to  assist  in  determining  the 
robustness  of  the  technique  with  decreasing  separation  distance  from  the  acoustic  source 
and  varying  receive  array  element  spacings  and  numbers.  This  technique  will  require  that 
the  receive  array  be  placed  outside  the  boundaries  imposed  by  the  enclosing  CHIEF 
sphere. 

D,  RECEIVE  ARRAY  MEASUREMENTS/SIMULATION 

This  step  involves  the  acquisition  (or  simulation)  of  complex  receive  response  at 
each  of  the  receive  array  elements.  This  step  requires  that  data  be  collected  for  each  of 
the  receive  array  elements  at  acoustic  device  rotation  angle  intervals  corresponding  to 
(360  degrees  /symmetry  order  of  the  enclosing  sphere).  The  sampling  (or  later 
interpolation)  of  this  data  at  the  prescribed  angles  allows  a  direct  relationship  between  the 
data  collected  and  the  defined  “slices”  of  the  enclosing  sphere  (established  by  the 
selected  symmetry  order). 

E.  COMPUTE  ORIGINAL  A  AND  B  MATRICES  VIA  CHIEF  FOR  MATLAB 

This  step  involves  the  computation  of  the  original  A  and  B  matrices  using 
“CHIEF  for  MATLAB”  tools  .  The  A  and  B  matrices  are  formed  from  the  surface 
integrals  in  CHIEF.  The  following  is  taken  from  documentation  provided  for  “CHIEF  for 
MATLAB”  by  S.E.  Forsythe  (NUWCDIVNPT):  “The  standard  problem  statement  for 
any  CHIEF  problem  is:  Ap  =  Bv  +  pinc,  where  A  is  the  matrix  formed  from  surface 
integrals  of  the  normal  derivative  of  pressure  over  the  surface  patches,  B  is  the  matrix 
formed  from  surface  integrals  of  the  pressure  over  the  surface  patches,  pinc  is  the  vector 
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of  incident  pressures  on  the  surface  patches  due  to  sources  in  the  region  of  interest,  p  is 
the  (to  be  solved  for)  vector  of  total  pressures  on  the  surface,  and  v  is  the  (to  be  solved 
for)  vector  of  velocities  on  the  surface.” 

CHIEF  for  MATLAB  has  the  capability  to  provide  the  A  and  B  matrices  and 
cease  further  computations  (to  solve  for  resulting  pressures  and  velocities).  This  is 
desired  for  our  far  field  extrapolation  technique  since  the  A  and  B  matrices  will  be 
augmented  with  other  factors  in  the  next  step.  The  primary  inputs  to  create  the  A  and  B 
matrices  depend  upon  the  acoustic  parameters  of  sound  speed,  density  of  the  test  medium 
and  frequency  of  interest  and  upon  the  geometric  configuration  of  the  surface  of  interest. 
In  our  case,  the  surface  of  interest  is  the  enclosing  sphere  from  the  first  step.  Thus,  given 
a  defined  enclosing  sphere,  the  matrices  A  and  B  can  be  formed. 

F.  COMPUTE  AUGMENTED  MATRICES  Aa  ,  Ba  AND  Ma 

This  next  step  begins  the  departure  from  the  usual  manner  in  which  the  CHIEF 
program  is  utilized.  That  is,  typical  CHIEF  applications  require  the  definition  of  a 
transducer  geometry  followed  by  the  characterization  of  a  pressure,  veloeity  or 
impedance  across  all  defined  patches.  Specification  of  one  of  these  parameters  (pressure, 
veloeity  or  impedance)  is  required  to  make  a  reduced  system  (from  the  original  defining 
equations  of  Ap  =  Bv  +  pinc  that  is  uniquely  solvable. 

In  this  thesis  methodology,  the  only  completely  defined  geometries  are 
represented  by  the  CHIEF  enclosing  sphere  (an  artificial  surface)  and  the  physical 
location  of  the  receive  array  elements;  the  geometry  of  the  device  under  test  is  not 
required  to  be  known  exactly  (with  the  exception  that  its  largest  dimension  must  be 
known).  This  flexibility  to  work  with  a  minimal  set  of  known  transducer  geometries  was 
a  desired  feature  of  the  methodology. 
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The  other  signifieant  differenee  in  this  CHIEF  usage  is  that  the  pressure  and 
veloeity  terms  at  eaeh  pateh  location  on  the  defined  CHIEF  enclosing  sphere  are  both 
unknown  values  that  must  be  solved  for  using  additional  information  constraints. 

A  least  squares  solution  is  the  ultimate  goal  of  the  methodology  and  relies  upon 
the  use  of  field  values  of  pressure  as  the  starting  point.  The  relation  between  the  CHIEF 
A  and  B  matrices  and  pressure  and  velocity  terms  form  the  fundamental  relationship  of 
0  =  -A  p  +  B  V  and  is  examined  first  below.  Since  this  system  by  itself  is 
underdetermined,  the  system  is  augmented  by  other  relationships  to  achieve 
overdetermined  status.  The  augmentation  relies  upon  the  CHIEF  A  and  B  matrices  and 
their  relationship  to  field  pressures;  reliance  is  also  made  upon  expected  adjacency 
conditions  of  pressure  and  velocity  between  patches. 

The  A  and  B  matrices  formed  in  the  last  step  provide  the  starting  point  to  develop 
the  following  relationship: 

Ma  =  -AaP  +  BaV 

where: 

Aa  =  Augmented  A  Matrix 
Ba  =  Augmented  B  Matrix 

p  =  pressures  at  all  “patches”  on  the  enclosing  sphere  (a  vector) 

V  =  velocities  at  all  “patches”  on  the  enclosing  sphere  (a  vector) 

Ma  =  Augmented  Field  Points  Matrix  (complex  responses  at  defined  points  in  the  test 
medium) 

The  augmentation  of  the  matrices  in  this  step  puts  the  measured  data  and  defining 
physical  equations  in  a  form  that  will  be  suited  for  a  later  least  squares  solution 
determination.  The  augmentation  of  the  A  and  B  matrices  is  examined  first  and  is 
described  as  follows: 
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1)  0  =  -A  p  +  B  V 


This  relationship  states  that  the  eomputed  A  and  B  matriees  (from  the  geometrieal 
and  aeoustie  eonsiderations  of  the  previous  step)  when  multiplied  by  their  respeetive 
pressures  and  veloeities  result  in  zero  net  pressure  on  the  surfaee  when  there  is  no  souree 
in  the  region  of  interest  as  required  by  the  Helmholtz  theorem. 

2)  pressure  at  a  field  point  =  delG  p  +  G  v 

This  defining  equation  relates  the  Greens  funetion  and  its  gradient  along  with 
veloeity  and  pressure  values  at  eaeh  physieal  location  on  the  enclosing  sphere  to  complex 
response  values  measured  at  (near)  field  test  points.  The  Green’s  function  and  its 
gradient  are  determined  via  geometrical  considerations  (ie:  the  physical  location  of  the 
field  points  with  respect  to  the  enclosing  sphere).  The  pressures  at  field  points  are  as 
measured  (or  simulated)  by  the  receive  array.  As  in  the  first  step  of  this  augmentation, 
this  defining  equation  utilizes  the  (to  be  solved  for)  pressures,  p,  and  velocities,  v,  along 
with  factors  that  can  be  determined  (  G  and  delG). 

3)  pressure/velocity  equivalence  at  adjacent  patches 

This  relationship  implies  that  due  to  the  relatively  small  size  of  the  patches 
comprising  the  surface  of  the  enclosing  sphere,  it  is  reasonable  to  assume  that  adjacent 
patches  do  not  differ  significantly  in  pressure  or  velocity  from  one  another.  They  can,  in 
fact,  be  considered  to  be  approximately  equal.  This  allows  a  physical  attribute  to  be 
associated  with  the  (to  be  solved  for)  pressures  and  velocities.  Enforcing  this  condition 
with  a  small  weight  in  the  least  squares  solution  prevents  the  extreme  case  where  non- 
realizable  pressure  and  velocity  swings  are  required  for  mathematical  realization 
Thus,  this  step  essentially  states  that  pi  =  pi+i  and  Vi  =  Vi+i,  where  i  represents  some  patch 
location  and  i+1  represents  the  patch  next  to  it  in  the  vertical  direction. 
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The  resulting  augmentation  of  the  A  and  B  matriees,  therefore,  yields  the 
following  constituents: 


Aa  Matrix 

Ba  Matrix 

-Ap 

Bv 

delG  p 

Gv 

Pi  =  Pi+i 

Vi  =  Vi+I 

Note  that  pressure  terms  are  related  to  the  Aa  matrix  while  velocity  terms  are 
related  to  the  Ba  matrix.  (The  details  of  the  actual  construction  protocols  of  these 
matrices  are  examined  later.) 

The  construction  of  the  augmented  field  point  array,  Ma,  is  now  examined.  The 
Ma  matrix  represents  the  expected  outputs  or  measured  values  corresponding  to  the 
constituents  of  the  augmented  Aa  and  Ba  matrices.  For  the  three  defined  parameters 
discussed  above,  Ma  takes  on  the  following  values: 

1) Ma=  0 

for  the  (-Ap+Bv)  terms 

2)  Ma  =  complex  response  measured  (simulated)  at  each  element  of  the  receive 

Array  for  the  (delG  p  +  G  v)  terms 

3)  Ma  =  0 

For  the  adjacent  equivalence  terms  of  pressure  and  velocity  represented  by 
(pi  -  pi+i)  and  (vi  -  Vi+i) 
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Since  a  least  squares  solution  teehnique  will  be  utilized  in  the  next  step,  an 
advantage  in  optimizing  the  solutionis  gained  by  scaling  the  matrix  constituents.  These 
re-scaled  matrix  eomponents  still  aecurately  represent  the  conditions  and  equalities  of  the 
set-up.  They  do,  however,  de-emphasize  the  wide  variation  in  magnitudes  of  the  various 
quantities  and  allow  the  solution  criteria  to  proeeed  based  upon  more  equal  weightings. 
This  re-sealing  is  also  appropriate  to  identify  as  putting  the  defining  equations  in 
dimensionless  form. 

For  example,  the  Green’s  Function  integral,  G,  for  a  partieular  set-up  returns  an 
absolute  maximum  value  of  573  while  the  gradient  returns  an  absolute  maximum  value  of 
3.82e-4.  Adjusting  the  Green’s  Function  integral  elements  by  first  dividing  all  G 
elements  by  rhoc  (1500000)  and  then  dividing  all  matrix  elements  of  (the  new)  G  and 
delG  by  3.82e-4  gives  matrix  elements  on  both  G  and  delG  that  have  maximum  values  of 
approximately  one. 

While  the  above  protocol  establishes  more  or  less  equal  weightings,  the 
methodology  does  allow  the  user  the  opportunity  to  place  specific  emphasis  on  the 
contributions  of  the  theoretieally  defined  Helmholtz  integral  relations  or  of  the 
empirieally  aequired  near  field  measurement  data.  This  is  embodied  in  the  alpha  and  beta 
terms  defined  below. 

For  this  effort,  the  comprising  matrices  are  scaled  by  the  following  faetors: 


Matrix 

Scaling  Factor 

A 

alpha 

B 

alpha/rhoc 

delG 

beta 

G 

beta/rhoe 

M 

beta 
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with: 

alpha  =  alphaO  /  (maximum  absolute  value  in  A) 
beta  =  betaO  /  (maximum  absolute  value  in  delG) 


and: 

alphaO  and  betaO  are  as  defined  by  the  user  to  emphasize  the  theoretieal  terms 
(Helmholtz  integrals  relations)  or  emphasize  the  empirieal  data  (near  field 
measurements),  respeetively,  in  the  final  solution. 

This  eompletes  the  augmentation  of  the  defining  matriees.  It  is  signifieant  to  note 
at  this  point  the  expeeted  size  of  the  augmented  matrix  eonstituents.  This  information  is 
tallied  below: 

Size  (rows  x  eolumns) 
numv  X  (sph  numv  *  sph  sym) 

numh  X  (sph_numv  *  sph  sym) 

(  sph  numv-l)  x  (sph  numv  *  sph  sym) 

(  sph_numv  +  numh  +  2*sph_numv  -2)  x  1 


Matrix 

A 

B 

delG 

G 

Pressure  equivalenee 
Veloeity  equivalenee 


The  matrix  sizes  presented  here  do  not  eonstitute  the  full  matrix,  rather  they  are 
representative  of  one  “sliee”  of  the  full  eirculant  matrix  .  The  full  matrix  size  is  larger 
by  a  faetor  of  sph  sym  for  eaeh  of  the  above  eonstituents. 
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The  full  Aa  and  Ba  matrices  are,  therefore,  of  size: 


(sph_sym  *  (  numh  +  3*sph_numv  -2)  )  x  (sph_numv  *  sph_sym) 

while  the  full  Ma  matrix  is  of  size: 

(  sph_sym  *  (  numh  +  3*sph_numv  -2)  )  x  1 

G.  SOLVE  FOR  PRESSURES  AND  VELOCITIES  VIA  LEAST  SQUARES 

This  step  uses  the  least  squares  solution  technique  to  compute  the  pressure  and 
velocity  values  at  each  of  the  enclosing  sphere  patches.  The  output  for  this  step  consists 
of  (numv  *  sph  sym)  pressure  values  and  the  same  number  of  velocities. 

The  least  squares  technique  allows  an  overdetermined  system  (“m”  equations  in 
“n”  unknowns,  with  m>n)  to  be  evaluated  in  an  optimized  manner.  This  optimized 
manner  results  in  values  of  x  that  minimize  the  quantity  |Ax  -  b|^  for  the  originating 
equation  of  Ax  =  b. 

The  vector  x  that  minimizes  |Ax  -  b|^  is  the  solution  to  the  normal  equations: 

A^  Ax  =  A^  b 

This  vector  x  =  (A^  A)‘^  A^  b  is  the  least  squares  solution  to  Ax  =b. 

The  implementation  of  the  least  squares  technique  for  this  thesis  effort  substitutes 
a  matrix  that  will  be  designated  as  the  AB_Matrix  for  A  and  the  full  Ma  matrix  for  b. 

The  AB_matrix  consists  of  the  full  matrices  Aa  and  Ba  that  have  been  interleaved  as 
follows: 

Aai  Bai  Aa2  Ba2  Aas  Ba3. . .  Aan  Ban;  with  n  =  (sph_numv  *  sph_sym) 
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The  AB_Matrix  is  actually  a  circulant  matrix  of  size: 

(sph_sym  *  (  numh  +  3*sph_numv  -2)  )  x  (2  *  sph_numv  *  sph_sym) 

The  justification  for  the  interleaving  of  the  Aa  and  Ba  components  will  be  examined  in 
detail  later. 

The  least  squares  solution,  x,  given  both  AB_Matrix  and  Ma  ,  is  a  vector  of 
alternating  pressures  and  velocities  (of  total  length:  2  *  sph_numv  *  sph  sym).  These 
values  correspond  to  the  pressures  and  velocities  at  each  patch  on  the  enclosing  sphere. 

H.  COMPUTE  Ma-new  MATRIX 

This  step  takes  the  pressures  and  velocities  computed  in  the  last  step  and 
substitutes  them  into  the  following  equation  to  compute  the  new  augmented  field  point 
matrix,  Ma-new  • 

M^a-new  Aa  p  +  Ba  V 

The  Aa  and  Ba  matrices  are  the  same  as  that  defined  previously  in  step  #5.  The  p 
and  the  v  represent  the  pressure  and  velocities  computed  via  the  last  step. 

This  newly  computed  matrix,  Ma-new  ,  can  then  be  used  to  evaluate  the  validity  of 
the  enclosing  sphere  pressures  and  velocities,  p  and  v,  and  whether  the  solution  is  likely 
to  be  properly  optimized.  This  is  accomplished  by  comparing  Ma-new  quantitatively  to 
the  original  Ma  matrix. 

Since  Ma-new  Utilized  a  data  set  derived  from  a  least  squares  solution  (the  p  and  v 
values),  it  is  not  expected  that  it  will  match  up  exactly  one  for  one  with  the  originating 
Ma  matrix.  There  are  important  considerations  to  be  made,  however. 
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Before  listing  this  criteria,  a  short  summary  of  the  constituents  of  the  Ma  matrix  is 
first  reiterated.  The  Ma  matrix  is  divided  up  into  four  sections. 

1)  The  first  section  has  values  corresponding  to  the  (-Ap  +  Bv)  calculations. 
Ideally  these  values  should  be  zero. 

2)  The  second  section  contains  the  complex  response  data  at  each  element  of  the 
receive  array  (these  are  near  field  data  points).  These  are  the  terms  of 

(delG  p  +  G  v). 

3)  The  third  section  contains  the  pressure  adjacency  terms  (pi  -  pi+i).  These 
values  should  be  close  to  zero. 

4)  The  fourth  section  contains  the  velocity  adjacency  terms  (vi  -  Vi+i).  These 
values  should  be  close  to  zero. 

The  following  are  the  primary  criteria  that  impact  the  successful  use  of  this 
technique  for  far  field  extrapolation. 

1)  The  complex  response  data  (section  two)  of  the  two  matrices,  Ma-new  and 
Ma,  should  agree  with  one  another  as  closely  as  possible. 

2)  The  (-Ap  +  Bv)  term  (section  one)  of  the  Ma-new  matrix  should  be 
minimized  as  much  as  possible  since  this  represents  a  known  physical  parameter. 

3)  The  pressure  adjacency  and  velocity  adjacency  terms  (sections  three  and 
four)  are  expected  to  be  small  terms,  but  this  condition  does  not  require  a  heavy 
enforcement.  Some  differences  in  pressure  and  velocity  values  on  adjacent 
patches  are  expected.  This  criteria  was  originally  added  to  lessen  the  chance  of 
skewing  the  pressure  and  velocity  values  to  physically  unrealizable  states. 
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4)  Since  “zeros”  are  not  likely  to  be  obtained  for  seetions  1,3  or  4,  it  is 
important  to  eonsider  the  difference  (in  dB)  between  the  magnitudes  in  these 
seetions  and  the  magnitudes  of  the  eomplex  response  data  at  each  element  of  the 
receive  array  (seetion  2).  This  differenee  should  be  maximized  as  mueh  as 
possible. 

The  data  presented  below  shows  the  two  matrices,  Ma-new  and  Ma,  for  a  given  set¬ 
up.  This  data  presented  has  been  optimized  to  a  reasonable  level  (as  examined  later)  and 
would  serve  to  qualify  the  solved-for  pressure  and  velocities  for  the  last  step  of  this 
methodology,  the  step  when  the  far  field  computations  are  aetually  made. 

The  four  distinct  sections  are  most  easily  observed  in  the  Ma-new  figure.  The 
sections  1  thru  4,  as  defined  above,  appear  in  succession  from  the  top  of  the  figure  to  the 
bottom.  The  figure  shows  magnitudes  expressed  in  dB;  the  colorbars  to  the  right  of  the 
figures  designate  the  correspondence  between  color  and  magnitude.  The  reddish  band  of 
colors  correspond  to  the  complex  response  data  (section  2)  and,  as  expected,  are  the 
highest  magnitudes  in  the  matrix. 


Figure  4  -  Original  and  “pv  solution”  Augmented  Field  Point  Matrices  (Ma  and  Ma-new) 
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I. 


COMPUTE  FIELD  POINT  RESPONSES 


This  last  step  uses  the  pressure  and  veloeity  data  obtained  and  evaluated  from  the 
previous  two  steps  to  compute  complex  response  at  defined  field  points.  This 
computation  proceeds  via  the  use  of  the  Greens  function  integral  and  its  gradient  as  has 
been  accomplished  in  earlier  steps  when  defining  the  augmented  matrices.  This  familiar 
relationship  is: 

p  =  (delG  p  +  G  v) 

As  before,  the  p  and  v  values  represent  the  pressures  and  velocities  at  each  patch 
location  on  the  enclosing  sphere.  The  G  and  delG  terms  represent  the  Greens  function 
integral  and  its  gradient  obtained  at  specified  field  points  with  respect  to  the  enclosing 
sphere.  The  selection  of  field  points  will  typically  be  made  to  encircle  the  acoustic 
device  on  a  designated  plane  allowing  what  are  traditionally  called  XY,  XZ,  or  YZ 
directivities. 

It  is  significant  to  note  here  that  there  is  no  distinction  in  this  step  whether  defined 
field  points  are  in  the  near  field  or  in  the  far  field.  This  technique  takes  a 
pressure/velocity  defined  enclosing  sphere  and  combines  its  elemental  traits  via  known 
physical  relationships  (Greens  function  integral  and  its  gradient)  to  achieve  pressure 
values  at  defined  locations.  Normalizing  directivity  data  obtained  via  this  technique  for 
an  increasing  range  will  show  the  expected  convergence  to  a  single  far  field  pattern  of 
defined  character. 
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V.  SOLUTION  SET  USING  METHODOLOGY 


This  section  takes  the  theory  and  methodology  detailed  in  the  previous  sections  to 
solve  a  real  world  problem. 

A.  DEFINITION  OF  SET-UP 

In  this  case  (and  as  supported  by  the  referenced  figures,  Table  1  and  the 
MATLAB  coding  provided),  a  horizontal  transmitting  line  array  of  length  1.885  meters  is 
operated  at  a  frequency  of  approximately  4775  Hz.  The  expected  near  field/far  field 
boundary  for  this  device  occurs  at  a  distance  of  approximately  11.32  meters.  This 
distance  exceeds  the  physical  size  of  the  Acoustic  Pressure  Tank  Facility  (APTF),  11.1 
meters  in  length,  where  this  device’s  performance  under  temperature  and  pressure 
conditions  are  desired  to  be  characterized.  It  is  an  ideal  candidate  to  investigate  the 
methodology  detailed  within  this  thesis. 

Figures  8  and  9  show  the  expected  theoretical  near  field  and  far  field  responses, 
respectively,  of  the  defined  device.  Figure  8  shows  near  field  performance  at  a  range  of 
approximately  5.341  meters.  This  test  distance  is  within  the  available  range  of  separation 
distances  that  can  be  achieved  within  the  APTF.  The  theoretical  far  field  response  was 
computed  at  a  test  distance  of  about  1131  meters.  This  places  device  performance 
characterization  well  into  far  field  conditions.  (Please  note  that  the  designation  of 
element  #6  on  the  near  field  data  ,  Figure  8,  corresponds  to  response  in  the  XY  plane.  As 
will  be  described  later,  element  #6  is  the  mid-point  element  of  the  1 1  element  vertical  line 
array  used  in  this  effort.) 

A  vertical  near  field  array  of  approximate  length  1.571  meters  was  used  for  this 
set-up.  This  near  field  array  is  made  up  of  1 1  equally  spaced  elements  and  was  placed  at 
a  distance  of  5.341  meters  from  the  origin.  (The  horizontal  transmitting  array  will  be 
rotated  about  this  same  origin  in  the  XY  plane  for  this  set-up.  This  origin  location 
corresponds  to  the  mid-point  of  the  transmitting  array.)  The  lengths  specified  for  both 
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the  near  field  array  and  for  the  transmitting  array  fit  within  the  confines  of  the  3.81 
diameter  APTF  test  vessel. 

The  selection  of  11  elements  for  the  near  field  array  was  arrived  a  bit  by 
experimentation  and  by  a  concerted  effort  to  reduce  the  numbers  of  acquisition  channels 
that  ultimately  would  be  required  for  this  implementation  (ie  -  16  channel  acquisition 
cards  are  fairly  easy  to  procure  as  off-the-shelf  components  and  would  be  suited  to  a 
future  hardware  implementation  for  this  technique). 

Figure  7  shows  a  section  of  the  enclosing  CHIEF  sphere  that  was  used.  The 
enclosing  sphere  was  configured  to  have  a  radius  of  about  1.384  meters.  The  sphere 
section  displayed  in  Figure  7  is  made  up  of  84  individual  patches.  A  symmetry  order  of 
180  was  selected;  this  yields  a  total  of  15,120  patches  on  the  enclosing  sphere  surface.  It 
is  a  part  of  this  methodology  to  solve  for  a  single  pressure  and  a  single  velocity  at  each 
one  of  these  patches.  Thus,  30,240  “unknowns”  must  be  solved  for  in  this  case.  (As  will 
be  discussed  later,  the  number  of  individual  patches  required  on  the  enclosing  surface  is 
the  primary  limiting  factor  for  this  methodology.  High  frequency  devices  with  large 
apertures  will  likely  require  an  enclosing  sphere  of  relatively  large  size  yielding  an 
extreme  number  of  patches.  This  likely  prohibits  solution  via  this  implementation  on 
standard  personal  computing  hardware  and  available  COTS  software.) 

B,  INITIAL  TEST  RESULTS 

Figures  10  through  14,  inclusive,  show  outputs  obtained  via  the  second  step  of  the 
MATLAB  coding.  Figure  10  shows  a  comparison  in  absolute  magnitude  of  the  original 
data  set  (obtained  via  simulation  in  this  example  but  would  correspond  to  actual 
measured  data  in  a  real  world  implementation)  and  the  hydrophone  near  field  response 
obtained  via  the  technique.  These  latter  hydrophone  responses  were  obtained  using  the 
solved  for  pressure  and  velocities  from  the  technique.  Figure  10  actually  shows  the 
difference  between  these  two  data  sets.  All  differences  are  less  than  |0.015|  dB  showing 
an  excellent  match. 
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The  remaining  figures,  11  thru  14,  allow  the  operator  to  assess  the  other 
eonditions  of  the  methodology.  Specitioally,  the  values  eomputed  for  the  Helmholtz 
Integral  relations  and  those  corresponding  to  the  pressure  and  velocity  adjacency 
conditions  need  to  be  assessed  as  well  to  judge  acceptance  of  the  solved  for  pressures  and 
velocities.  In  an  ideal  solution,  these  other  conditions  would  be  equal  to  zero.  Since  the 
least  squares  solution  technique  was  utilized,  a  best  fit  to  the  defining  equations  was 
determined  based  upon  near  field  array  element  data  and  assumed  zero  conditions  for  the 
Helmholtz  Integral  relations  and  pressure/velocity  adjacency  conditions.  Is  it  not  a  true 
requirement  to  have  these  latter  two  conditions  exactly  equal  to  zero.  It  is,  however, 
desirable  to  have  these  other  conditions  much  smaller  than  the  computed  hydrophone 
responses.  It  is  for  this  reason  that  dB  levels  are  reported  in  Figures  1 1  through  14. 

Figure  1 1  shows  a  plot  of  the  minimum  hydrophone  responses,  maximum  ZeroO 
and  maximum  Zerol  conditions  on  a  per  angle  basis.  The  nomenclature  ZeroO  and  Zerol 
correspond  to  the  “expected  zeros”  for  the  Helmholtz  Integral  relations  and 
pressure/velocity  adjacency  conditions,  respectively.  Since  the  need  here  is  to  assess  the 
solved  for  pressures  and  velocities,  it  is  desirable  to  look  at  worst  case  conditions  in  these 
data  sets.  That  is,  the  selection  of  minimum  hydrophone  response  and  maximum  ZeroO 
and  Zerol  responses  provides  the  smallest  difference  between  the  two  data  set  conditions. 
(Reference  Figure  3  -  Original  and  “pv  solution”  Augmented  Field  Point  Matrices  (Ma 
andMa-new)  for  visualization  of  the  ZeroO  and  Zerol  conditions.  ZeroO  values  are  shown 
above  the  reddish  band,  the  complex  response  data,  while  Zerol  values  are  shown 
below.) 


Figures  12  and  13  augment  Figure  11  by  providing  the  actual  computed 
differences  described  above.  In  this  data  set,  a  minimum  difference  on  the  order  of  about 
44  dB  and  54  dB  are  observed  for  the  ZeroO  and  Zerol  conditions.  These  magnitudes, 
from  experimentation,  are  likely  to  indicate  that  the  solved  for  pressures  and  velocities 
are  acceptable  to  use  for  far  field  extrapolation.  The  adjustment  of  the  values  for  the 
alphaO  and  betaO  terms  will  be  looked  at  later  to  try  and  better  this  solution. 
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Figure  14  is  perhaps  the  most  useful  of  all  of  the  figures  in  the  Figure  1 1  through 
14  set.  Figure  14  provides  a  “snapshof’  of  all  of  the  parameters  discussed  above 
(Helmholtz  Integral  Relations,  Near  Field  Hydrophone  Response,  Pressure  Adjacency 
and  Velocity  Adjacency  conditions).  It  provides  a  means  to  directly  compare  these 
parameters  with  their  originating  values.  It  also  allows  the  operator  to  make  direct 
comparison  amongst  the  Ma-new  constituents  (which  are  the  parameters  mentioned  above). 
The  color  in  this  figure  is  representative  of  magnitude  as  shown  in  the  colorbar  scale 
provided.  Thus,  it  is  possible  to  roughly  deduce  at  a  glance,  all  of  the  differences  as  noted 
in  Figure  1 1  through  13. 

Figures  15  through  18  used  the  solved  for  pressures  and  velocities  obtained  for 
alphaO=l  and  betaO=l  (with  the  performance  characteristics  as  identified  in  the  previous 
figures)  to  solve  for  pressure  at  defined  field  points.  These  defined  points  correspond  to 
the  original  near  field  points  and  points  in  the  far  field  (which  is  the  motivation  for  this 
technique).  Figures  15  and  16  show  that  the  solved  for  pressures  and  velocities  yield  near 
field  pressures  that  have  an  extremely  high  degree  of  agreement  with  the  original  near 
field  data  set.  Absolute  differences  of  less  than  .008  dB  were  observed. 

Figures  17  and  18  show  far  field  response  as  originally  simulated  and  as 
extrapolated  by  this  technique.  The  majority  of  the  differences  in  these  two  data  sets 
(Figure  18)  are  on  the  order  of  a  few  tenths  of  a  dB  with  just  a  few  outliers  (these  are 
addressed  later).  Thus,  the  technique  has  done  an  excellent  job  of  matching  expected  far 
field  performance  with  that  obtained  via  the  use  of  solved  for  pressures  and  velocities  on 
an  enclosing  sphere.  The  technique  is  therefore  viable. 

As  mentioned  above,  the  majority  of  the  differences  in  the  two  data  sets  (Figure 
18)  are  on  the  order  of  a  few  tenths  of  a  dB.  There  are,  however,  some  difference  values 
as  high  as  about  4.3  dB.  Closer  examination  of  the  locations  of  these  large  differences 
shows  a  correspondence  with  a  null  in  the  directional  response.  Thus,  the  technique  did 
an  excellent  job  of  producing  a  far  field  response  but  missed  representing  to  a  high 
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accuracy  the  depth  of  the  nulls  in  these  far  field  responses.  This  will  be  looked  at  again 
onee  new  values  of  alphaO  and  betaO  are  utilized. 

The  values  for  alphaO  and  betaO  are  now  ehanged  to  evaluate  their  effeet  on  the 
solved  for  pressures  and  veloeities  and  eventually  their  effeet  on  extrapolated  far  field 
performanee. 

C.  SENSITIVITY  OF  SOLUTION  TO  VARIATION  IN  BETAO 

Figures  19  through  23  used  values  of  alphaO=l  and  beta0=10.  (Reeall  that  the 
betaO  term  emphasizes  near  field  element  data.)  Figure  23  shows  some  improvement  in 
the  differenee  in  level  between  the  hydrophone  response  measured  and  its  other  zeroO 
and  zerol  eonstituents.  These  minimum  differenees  are  now  about  74  and  64  dB, 
respeetively.  Pressure  eomputations  at  the  various  field  points  were  not  eomputed. 

Figures  24  through  28  used  values  of  alphaO=l  and  beta0=100.  An  improvement 
in  differenee  levels  was  noted  (97  dB  and  84  dB) ,  however,  pressure  eomputations  were 
not  undertaken. 

Figures  29  through  33  used  values  of  alphaO=l  and  beta0=500.  Differences  of 
110  dB  and  98  dB  between  hydrophone  response  and  ZeroO  and  Zerol  eonditions, 
respeetively,  were  noted.  Figure  33  shows  this  same  level  of  differenee  by  inspeetion. 

Despite  ehanges  in  magnitude  between  these  various  Ma-new  figures,  an 
examination  of  Figure  33  and  Figure  14  shows  how  the  betaO  term  is  foreing  a  larger 
differenee  between  the  hydrophone  response  data  and  the  Helmholtz  Integral  relations 
and  pressure/veloeity  adjacency  conditions.  Each  of  these  Ma-new  figures  has  a  dynamic 
range  of  150  dB.  Figure  33  shows  how  the  hydrophone  response  data  is  being 
emphasized  while  the  other  parameters  are  being  suppressed. 

Computation  of  pressure  at  the  various  field  points  was  undertaken  for  the 
alpha0=l  and  beta0=500  eonditions.  Figures  34  through  37  show  the  resulting  pressure 


33 


responses  and  simulated/teehnique  derived  differenees  for  these  values.  As  before,  the 
near  field  data  eontinues  to  show  solid  agreement  between  simulated  and  teehnique 
derived  values  (Figures  34  and  35).  Figures  36  and  37  show  far  field  performanee  traits. 
The  noteworthy  part  of  this  data  set  is  that  the  differenees  between  the  two  far  field  data 
sets  are  now  less  than  1  dB  (down  from  4.3  dB  in  Figure  18).  Thus,  this  betaO  value 
ehange  has  resulted  in  pressures  and  veloeities  that  provide  a  better  far  field  fit.  This  is 
espeeially  seen  now  in  the  better  agreement  afforded  at  the  nulls  of  the  far  field  response. 

A  final  betaO  ehange  to  betaO  =  1000  was  undertaken.  Figures  38  through  43 
show  data  obtained  with  alpha0=l  and  beta0=1000.  As  ean  be  seen  on  the  far  field 
differenee  reported  in  Figure  43,  there  was  little  ehange  in  far  field  response  provided  by 
this. 


While  far  field  agreement  in  this  ease  eould  have  essentially  stopped  with 
alpha0=l  and  beta0=l,  the  other  values  of  betaO  were  examined  to  evaluate  the 
teehnique.  Inereases  in  betaO  did  lead  to  a  better  fit  in  the  nulls  of  the  far  field,  but  this  is 
likely  of  limited  praetieal  value. 

D,  SENSITIVITY  OF  SOLUTION  TO  VARIATION  IN  ALPHAO 

In  likewise  fashion,  to  further  evaluate  the  teehnique,  ehanges  in  the  value  of 
alphaO  (emphasizing  the  theoretieal  portion  of  the  teehnique)  will  now  be  undertaken. 

Figures  44  through  48  show  data  eorresponding  to  alpha0=10  and  beta0=l. 

Figures  49  through  53  show  data  eorresponding  to  alpha0=100  and  beta0=l. 

Figures  54  through  58  show  data  eorresponding  to  alpha0=100  and  beta0=l. 

An  examination  of  the  Ma-new  figures  for  eaeh  of  these  eonditions  shows  that  there 
is  little  to  no  ehange  in  level  differenees  between  hydrophone  response  and 
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pressure/velocity  adjacency  conditions.  The  Helmholtz  Integral  relations  are  showing  a 
suppression  as  alphaO  increases. 

Figures  59  through  62  show  pressure  responses  computed  for  the  near  field  and 
far  field  cases.  An  examination  of  the  far  field  difference  figure  (Figure  62)  shows  that 
very  little  was  gained  by  increasing  the  alphaO  term.  The  differences  in  the  nulls  that 
were  noted  for  the  alphaO=l  and  betaO=l  far  field  computations  still  exist  and  are  at  the 
same  approximate  level. 

E,  SUMMARY  OF  SENSITIVITY  OF  SOLUTION  TO  ALPHAO  AND  BETAO 

Table  1  summarizes  the  conditions  examined  in  the  previous  examples.  The 
following  definition  of  terms  applies  to  the  data  presented. 

alphaO: 

This  is  the  weighting  term  to  emphasize  the  theoretical  Helmholtz  Integral 
Relations. 

betaO: 


This  is  the  weighting  term  to  emphasize  empirical  data 

MAX  |Near  Field  Hydro  Diflj: 

This  value  represents  the  maximum  difference  between  the  measured  hydrophone 
responses  of  all  near  field  receive  array  elements  and  the  values  computed  by  the  second 
step  of  the  methodology.  A  low  value  represents  good  agreement  and,  therefore,  a 
possible  candidate  for  far  field  extrapolation  (using  the  pressure  and  velocity  values  so 
obtained). 


35 


MIN  min  hydro  minus  max  zeroO: 

This  value  represents  the  minimum  differenee  between  the  minimum  near  field 
hydrophone  data  and  the  maximum  zeroO  data  (theoretieal  Helmholtz  relations)  obtained 
at  eaeh  angular  position.  A  large  magnitude  here  also  signifies  that  the  eonditions  are 
favorable  for  far  field  extrapolation. 

MAX  min  hydro  minus  max  zeroO: 

This  value  represents  the  maximum  differenee  between  the  minimum  near  field 
hydrophone  data  and  the  maximum  zeroO  data  (theoretieal  Helmholtz  relations)  obtained 
at  eaeh  angular  position.  This  is  largely  for  informational  purposes  only. 

MIN  min  hydro  minus  max  zerol: 

This  value  represents  the  minimum  differenee  between  the  minimum  near  field 
hydrophone  data  and  the  maximum  zerol  data  (pressure  and  veloeity  equivaleney  for 
adjaeent  patehes)  obtained  at  eaeh  angular  position.  A  large  magnitude  here  also 
signifies  that  the  eonditions  are  favorable  for  far  field  extrapolation. 

MAX  min  hydro  minus  max  zerol: 

This  value  represents  the  maximum  differenee  between  the  minimum  near  field 
hydrophone  data  and  the  maximum  zerol  data  (theoretieal  Helmholtz  relations)  obtained 
at  eaeh  angular  position.  This  is  also  largely  for  informational  purposes  only. 

MAX  |Far  Field  Difference! : 

This  value  represents  the  maximum  differenee  between  the  methodology- 
eomputed  far  field  response  and  a  measured  (or  in  this  ease,  simulated)  far  field  response 
obtained  at  eaeh  angular  position.  A  low  value  indieates  that  the  methodology  produeed 
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a  reasonable  approximation  to  actual  far  field  performance.  A  high  value  may  indicate  an 
overall  poor  match  or  perhaps  just  a  poor  representation  at  certain  angles  (most  likely  in 
the  nulls).  It  is  advantageous  to  look  at  the  overall  pattern  obtained  to  further  define  the 
severity  of  a  high  value.  (Please  note,  however,  that  all  constituents  to  obtain  this  data 
value  would  not  typically  be  available  in  practice  (ie  -  the  purpose  of  this  thesis  was  to 
achieve  a  far  field  formulation  when  only  near  field  test  conditions  were  available). 
These  results  were  available  in  this  thesis  effort  and  therefore  have  been  examined.) 

The  initial  equal  weighting  of  the  alphaO  and  betaO  terms  (by  setting  both  equal  to 
1)  for  this  example  gave  excellent  results  with  little  gained  by  any  other  adjustments 
(other  than  a  bit  better  definition  of  the  response  in  far  field  nulls).  This  is  not  necessarily 
likely  to  be  accurate  in  all  cases.  Where  there  is  some  uncertainty  in  the  measured  data 
acquired,  a  better  fit  may  be  obtained  by  emphasizing  the  theoretical  portion  of  the 
technique.  This  will  be  examined  in  the  next  section  when  error  in  range  and  positioning 
of  the  near  field  array  are  examined. 
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MAX 


MIN 


MAX 


alphaO  betaO 


INear  Field  Hydro  Diff| 


min_hydro_minus_max_zeroO 


min_hydro_minus_max_zeroO 


(dB) 


(dB) 


(dB) 


1 

1 

0.010961319 

54 . 0 

1 

10 

0 . 001350528 

74.2 

1 

100 

0.000097393 

96.7 

1 

500 

0.000008493 

110 . 5 

1 

1000 

0.000009353 

116.5 

10 

1 

0.011399652 

71 . 7 

100 

1 

0 . 011405154 

91 . 7 

500 

1 

0 . 011405208 

105 . 7 

93.9 

113.9 
132 . 6 
146.2 
152 . 1 
111 .  4 

131 . 4 

145 . 4 


MIN  MAX  MAX 

alphaO  betaO  min_hydro_minus_max_zerol  min_hydro_minus_max_zerol  |  Far  Field  Difference! 


(dB) 

(dB) 

(dB) 

1 

1 

43.8 

89.2 

4.30 

1 

10 

63.8 

107 . 6 

n/a 

1 

100 

84.2 

125.2 

n/a 

1 

500 

97  .  9 

138.9 

0.95 

1 

1000 

103.8 

144 . 8 

0.90 

10 

1 

43.6 

88.7 

n/a 

100 

1 

43.6 

88.6 

n/a 

500 

1 

43.6 

88.6 

4.00 

Table  1  - 

Solution  Set  Methodology  Data  Summary 
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VI.  ERROR  ANALYSIS  -  EFFECT  OF  POSITIONAL  UNCERTAINTY 


A,  INTRODUCTION 

This  section  addresses  positional  uncertainties  in  the  placement  of  the  near  field 
array  and  its  resulting  effect  on  the  extrapolated  far  field  response  using  the  described 
methodology.  The  examples  provided  herein  were  computed  using  the  same  device 
configurations,  acoustic  parameters  and  CHIEF  set-ups  detailed  in  the  previous  Solution 
Set  using  Methodology  section.  In  this  section,  however,  the  alphaO  and  betaO  terms 
were  established  at  unity  each  and  were  not  allowed  to  vary. 

Positional  uncertainty  was  defined  to  correspond  to  an  inexact  determination  of 
range,  vertical  placement  or  angular  orientation  of  the  near  field  array  with  respect  to  the 
transmitting  array’s  origin.  The  data  for  this  section  was  obtained  by  varying  one  of  the 
aforementioned  parameters  and  then  running  through  the  computation  steps  in  the 
methodology  to  achieve  extrapolated  far  field  response  performance.  The  assessment  of 
the  error  so  induced  is  primarily  based  upon  the  examination  of  the  far  field  responses 
and  their  differences.  Near  field  responses  (showing  excellent  agreement  throughout  the 
various  error  conditions)  and  figures  of  the  Ma-new  matrices  are  provided  for 
completeness  in  this  section. 

B,  RANGE  ERROR 

Figures  63  through  86,  inclusive,  examine  the  effect  of  intentionally  induced 
errors  in  range.  For  these  figures,  the  simulated  data  used  range  values  with  error  while 
the  thesis  technique  used  the  nominal  range  magnitudes  (ie:  range  without  error). 
Various  range  errors  were  utilized  to  assist  in  evaluating  the  sensitivity  of  the  technique 
to  this  parameter.  These  are  summarized  on  the  following  page. 
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Figures  63  through  68 


-2*lambda  (-0.628  meters) 


Figures  69  through  74  +2*lambda  (+0.628  meters) 

Figures  75  through  80  -1  em 

Figures  81  through  86  +1  em 

For  referenee,  Figures  10  through  18  eorrespond  to  a  range  error  equal  to  zero 
(and  with  alpha0=l  and  beta0=l). 

As  expeeted,  the  larger  range  errors  provided  the  most  error  in  the  extrapolated  far 
field  response  as  shown  in  the  far  field  difference  data  reported  in  Figures  68  and  74. 
Even  in  these  extreme  cases  (+/-  0.628  meters  on  a  defined  range  of  about  5.354  meters), 
however,  the  far  field  response  (neglecting  the  response  in  the  nulls  of  the  response) 
differed  from  the  zero  range  error  case  on  average  no  more  than  about  0.5  dB.  These 
results  were  a  bit  unexpected. 

The  range  errors  of  +/-  1  cm  were  computed  in  this  data  set  to  better  represent 
range  errors  that  are  more  likely  to  be  encountered  within  the  test  facility  (APTF)  . 
These  error  conditions  (Figures  74  and  80)  show  very  little  difference  from  the  zero  range 
error  case.  Thus,  expected  uncertainties  in  the  range  (to  +/-  1  cm  and  for  the  set-up 
conditions  examined)  are  not  deemed  to  be  a  significant  source  of  error  in  extrapolating 
far  field  response  via  the  technique  presented. 

C.  VERTICAL  PLACEMENT  ERROR 

Figures  87  through  104,  inclusive,  examine  the  effect  of  intentionally  induced 
errors  in  the  vertical  placement  of  the  receive  array.  In  like  fashion  as  for  the  range  error 
case,  the  simulated  data  used  the  error  condition  while  the  thesis  technique  used  the 
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nominal  values  (ie;  without  error).  Various  errors  were  utilized  to  assist  in  evaluating  the 
sensitivity  of  the  technique  to  this  parameter.  These  are  summarized  below; 

Figures  87  through  92  +lambda/2  (+0.157  meters) 

Figures  93  through  98  +lambda  (+0.314  meters) 

Figures  99  through  104  -lambda/2  (-0.157  meters) 

The  selection  of  these  error  values  as  multiples  of  lambda/2  stems  from  the 
configuration  of  the  near  field  receive  array  itself.  In  the  example  that  has  been 
discussed,  the  spacing  of  the  elements  in  the  near  field  array  is  lambda/2.  The  worst  error 
condition  that  was  deemed  likely  to  occur  in  practice  was  considered  to  be  an  offset  in  the 
placement  of  the  array  by  one  element.  A  2  element  offset  (+lambda  in  this  example) 
was  selected  to  further  assess  the  technique.  This  latter  vertical  error  is,  however,  not 
likely  to  occur  in  practice  for  the  inter-element  spacings  specified  in  this  example.  As  for 
the  range  error  case.  Figures  10  through  18  can  be  referenced  for  comparison  to  a  zero 
error  case. 

Again,  as  expected,  the  larger  vertical  placement  errors  provided  the  most  error  in 
the  extrapolated  far  field  response  as  shown  in  the  difference  data  reported  in  Figure  98. 
For  the  more  typical  error  cases  (vertical  errors  of  +/-  lambda/2),  the  differences  in  far 
field  response  were  comparable  with  those  obtained  for  the  zero  error  case.  None  of  the 
vertical  placement  errors  accounted  for  more  than  perhaps  a  0.25  dB  error  in  the  final 
reported  far  field  response  (again,  neglecting  the  response  in  the  nulls  of  the  response). 

Thus,  uncertainties  in  the  vertical  placement  of  the  near  field  array  (to  within  an 
element  spacing  and  for  the  set-up  conditions  examined)  are  not  deemed  to  be  a 
significant  source  of  error  in  extrapolating  far  field  response  via  the  technique  presented. 
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D,  ANGULAR  ORIENTATION  ERROR 


Figures  105  through  110,  inclusive,  examine  the  effect  of  intentionally  induced 
errors  in  the  angular  orientation  (offset)  of  the  receive  array  with  respect  to  the 
transmitting  array.  In  like  fashion  as  for  the  previous  two  cases,  the  simulated  data  used 
the  error  condition  while  the  thesis  technique  used  the  nominal  values  (ie:  without  error). 
A  single  value  of  error  equivalent  to  a  10  degree  offset  was  examined  in  this  data  set. 

Figure  110  shows  the  difference  between  the  simulated  data  set  and  that  obtained 
from  the  thesis  technique.  Upon  initial  inspection,  this  figure  seems  to  indicate  gross 
errors  caused  by  the  10  degree  displacement.  Differences  of  close  to  30  dB  are  noted  in 
this  figure.  Upon  closer  examination  of  the  defining  constituents  of  the  difference  figure 
(the  data  presented  in  Figurel09),  it  is  readily  seen  that  the  shape  and  structure  of  the  two 
far  field  responses  are  essentially  the  same.  The  introduction  of  the  10  degree  angular 
error  has  shifted  the  response  accordingly  and  caused  the  large  difference  errors.  While 
this  results  in  large  difference  errors  (Figure  110),  this  condition  is  easily  post-processed 
and  does  not  represent  any  substantial  source  of  error  or  a  deterrent  to  the  use  of  the 
thesis  methodology. 

Thus,  as  for  the  other  cases,  uncertainties  in  the  angular  orientation  of  the  near 
field  array  are  not  deemed  to  be  a  significant  source  of  error  in  extrapolating  far  field 
response  via  the  technique  presented.  This  statement  is  likely  to  hold  true  for  all  set-ups 
and  not  just  for  the  unique  parameters  looked  at  in  this  example  set. 

E,  ERROR  SUMMARY  COMMENTS 

Table  2  summarizes  the  error  conditions  examined  in  the  previous  examples  in 
this  section.  Please  reference  the  previous  section.  Solution  Set  using  Methodology,  for 
definitions  of  the  nomenclatures  utilized. 
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For  the  given  example  set-up,  errors  in  range,  vertieal  plaeement  and  angular 
orientation  of  the  near  field  reeeive  array  (when  taken  singularly  and  not  taken  in  any 
kind  of  an  extreme  sense)  indueed  very  little  error  in  the  extrapolated  far  field  response 
(again,  negleeting  the  nulls  of  the  response)  Attention  to  error  quantifieation  in  terms  of 
some  lambda  value  (or  fraetional  value)  was  given  in  an  effort  to  assist  in  extrapolating 
this  performanee  data  to  other  set-up  eonditions.  While  this  ability  to  extrapolate 
performanee  traits  is  likely  aehievable,  the  ability  to  position  deviees  within  the  same 
lambda  toleranee  will  be  hampered  as  operating  frequeneies  inerease  (and  wavelengths 
eorrespondingly  deerease). 

The  overall  eonelusion  provided,  however,  is  that  the  teehnique  worked  extremely 
well  at  charaeterizing  far  field  response  even  when  physieal  set-up  parameters  were 
known  to  be  in  error. 
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alphaO=betaO=l 


MAX 


MIN 


MAX 


error  condition 
no  error 

range  error:  -2*lambda 
range  error:  +2*lambda 
range  error:  -1  cm 

range  error:  +1  cm 
vert  rec  array  error : 

+ lambda/ 2 

vert  rec  array  error:  ilambda 
vert  rec  array  error:  - 
lambda/ 2 

azimuthal  error:  +10  degrees 


error  condition 
no  error 

range  error:  -2*lambda 
range  error:  +2*lambda 
range  error:  -1  cm 

range  error:  +1  cm 
vert  rec  array  error: 

+ lambda/ 2 

vert  rec  array  error:  llambda 
vert  rec  array  error:  - 
lambda/ 2 

azimuthal  error:  +10  degrees 


I  Near  Field  Hydro  Diff| 
(dB) 

0.010961319 
0 . 015705251 
0 . 010182845 
0 . 011011699 
0 . 010911516 

0 . 019242789 
0 . 023584640 

0 . 019242789 
0 . 011221058 

MIN 

min_hydro_minus_max_zerol 

(dB) 

43.80 
40.43 

41 . 85 
43.83 
43 . 78 

35.09 

29.10 

35.09 

43.80 


min_hydro_minus_max_zeroO 

(dB) 

53.95 

49.88 
53.68 

53.96 
53.95 

50.60 

40 . 73 

50.60 
53 . 95 

MAX 

min_hydro_minus_max_zerol 

(dB) 

89.16 

83.27 

92 . 73 
89.05 
89.26 

87.88 
81.83 

87.88 

89.16 


min_hydro_minus_max_zeroO 

(dB) 

93.86 

88.29 

93.62 

93.78 

93.93 

89.97 

85.23 

89.97 

93.86 

MAX 

Far  Field  Difference! 
(dB) 

4.30 

12.60 
5.20 
4 .40 
4.25 

3.00 

2.60 

3.00 

30.00 


Table  2  -  Effect  of  Positional  Uncertainty  Data  Summary 
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VII.  ADDITIONAL  SOLUTION  SETS 
A,  INTRODUCTION 

This  section  addresses  the  effectiveness  of  the  technique  when  the  near  field  array 
is  placed  closer  to  the  device  under  test.  The  same  transmitting  array,  1.885  meters  in 
length  and  operating  at  a  frequency  of  about  4775  Hz,  is  used  in  all  cases.  The  cases 
presented  in  this  section  are  summarized  below: 


Near  Field  Array  Location 

alphaO 

betaO 

Figures  111  to  113 

2.082  meters 

1 

1000 

Figures  1 14  to  116 

1 .466  meters 

1 

1000 

Figures  117  to  119 

1 .466  meters 

1 

10,000 

*Figures  120  to  122 

1 .466  meters 

1 

1000 

*  -  Figures  120  to  122  represent  data  taken  with  a  different  configuration  near  field 
array.  The  array  used  for  this  data  set  was  about  2.199  meters  long  and  was  composed  of 
15  equally  spaced  elements. 

The  data  presented  here  was  provided  more  for  academic  reasons  than  for  any 
required  need.  The  location  of  the  near  field  array  at  a  range  of  5.354  meters  in  the 
previous  sections  adequately  met  the  restrictions  imposed  by  the  dimensions  of  the  APTF 
test  vessel.  Thus,  the  test  set-ups  specified  in  this  section  may  not  be  required,  but  they 
do  serve  to  address  the  capabilities  of  the  technique. 

B.  NEAR-FIELD  ARRAY  RANGE  OF  2,082  METERS 

Figures  111  through  113  used  a  near  field  array  range  of  2.082  meters.  The  near 
field  data  sets  (simulated  and  from  the  technique)  matched  very  well;  the  far  field  data 
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sets  also  matched  extremely  well.  It  is  interesting  to  note  that  the  near  field  data  actually 
now  shows  a  local  null  at  the  same  angular  location  where  the  far  field  data  actually 
peaks. 


Please  note  that  the  enclosing  sphere  required  for  these  data  sets  was  also  moved 
to  a  closer  proximity  to  the  device  under  test.  Since  the  radius  of  this  enclosing  sphere 
was  decreased  and  yet  the  frequency  of  operation  stayed  the  same,  fewer  numbers  of 
patches  were  required  in  the  vertical  direction  and  a  lower  symmetry  order  sufficed  for 
the  data  set.  The  lower  symmetry  order  (since  it  is  directly  related  to  the  number  of 
discrete  angles  in  the  responses)  accounts  for  the  choppiness  of  these  responses  when 
compared  to  the  previous  solution  sets’  figures. 

C.  NEAR-FIELD  ARRAY  RANGE  OF  1,466  METERS 

Figures  1 14  to  122  all  used  a  near  field  array  range  of  1.466  meters.  Figures  1 14 
through  116  and  Figures  117  through  119  used  the  same  alphaO  values  (one)  but  differing 
betaO  values  of  1000  and  10,000,  respectively.  The  far  field  response  of  Figure  115 
showed  a  very  reasonable  fit  to  the  simulated  data  along  the  main  beam  of  the  device,  but 
did  not  have  as  good  a  match  in  the  side  lobes.  The  data  for  Figure  118  was  computed 
with  a  higher  betaO  value  to  attempt  to  correct  this  side  lobe  problem.  The  Ma-new  data  for 
these  two  conditions  (Figures  116  and  119)  would  lead  one  to  believe  that  the  higher 
betaO  value  would  give  a  more  accurate  description  of  the  devices  actual  far  field 
response.  This  was  not  found  to  be  the  case.  The  far  field  response  for  these  two 
conditions  (Figures  115  and  118)  looked  to  be  very  comparable. 

After  some  experimentation,  it  was  observed  that  increasing  the  length  of  the  near 
field  array  yielded  an  improvement  in  side  lobe  performance  of  the  device  under  test. 
Figures  120  to  122  show  response  with  a  near  field  array  that  has  been  lengthened  by 
2*lambda  to  2.199  meters  (originally  this  array  was  about  1.571  meters  in  length).  The 
far  field  response  (Figure  121)  shows  an  improved  side  lobe  performance  over  the 
previous  data  sets  using  the  same  near  field  array  range  of  1 .466  meters. 
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Please  note  that  the  number  of  elements  in  the  2.199  meter  near  field  array  was 
inereased  from  11  to  15  in  an  effort  to  maintain  lambda/2  spaeing.  The  improvement  in 
far  field  response  observed,  however,  was  due  to  its  length  inerease  and  not  due  to  the 
inerease  in  the  number  of  its  comprising  elements.  This  latter  characteristic  was  arrived 
at  empirically  but  is  not  supported  with  any  data  in  this  write-up. 

D.  SUMMARY 

Based  upon  the  data  sets  presented  in  this  section,  it  is  desirable  to  quantify  a 
minimum  near  field  array  length  and  a  minimum  test  distance  for  a  given  test  set-up  to 
provide  some  assurance  that  proper  far  field  performance  is  extrapolated.  Using 
empirical  results,  it  is  believed  that  a  near  field  array  with  a  length  that  exceeds  the 
maximum  dimension  of  the  device  under  test  placed  at  a  separation  distance  comparable 
to  its  length  could  achieve  this  condition.  This  latter  hypothesis  is  speculated  but  is  not 
stated  as  proof  in  this  thesis.  It  is  however,  supported  by  the  data  presented  in  Figures 
123  through  125  (with  the  transmitting  array,  near  field  array  and  range  all  specified  at 
1.885meters). 
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VIII.  SUMMARY 


A,  INTRODUCTION 

A  technique  to  extrapolate  far  field  performance  of  an  acoustic  device  from  a  data 
set  obtained  in  the  near  field  was  demonstrated.  The  technique  successfully  married  the 
Helmholtz  integral  relationships  (A  and  B  matrices  and  pressure  computations  utilizing 
the  integration  of  the  Green’s  function  and  its  gradient)  and  expected  physical  parameters 
with  near  field  acquired  data  to  yield  an  over  determined  system.  Solution  via  a  least 
squares  technique  provided  pressure  and  normal  velocity  terms  for  each  patch  of  the 
enclosing  sphere  surface.  This  characterization  of  the  enclosing  sphere’s  surface  allowed 
extrapolation  to  any  range  (near  or  far  field)  when  Greens’  function  and  its  gradient 
integration  values  appropriate  to  that  range  are  utilized. 

B,  SOLUTION  SETS 

An  example  solution  set  was  provided  that  demonstrated  the  ability  of  the 
technique  to  produce  extremely  accurate  results.  Figures  4  and  5  below  show  the  data 
computed  for  the  near  field  and  far  field,  respectively,  for  the  defined  acoustic  device  (a 
1.885  meter  long  linear  array  operating  at  about  4775  Hz).  (Figures  4  and  5  represent  the 
same  data  set  shown  in  Figures  34  and  36.) 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20 ;  range  =  5,3541  meters;  array  length  =  1 .885  meters 
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Figure  5  -  Near  Field  Data  -  eomputed  and  simulated 


FAR  FIELD  DATA  (on  axis)  with  l<  =  20 ;  range  =  1130.9734  meters;  array  length  =  1.885  meters 


Figure  6  -  Far  Field  Data  -  eomputed  and  simulated 


The  example  set-up  was  speeifieally  designed  for  a  minimal  set  of  elements  in  the 
near  field  array  (eleven).  This  minimal  set  was  utilized  in  an  attempt  to  ease  the 
adaptation  of  the  technique  to  a  future  hardware  implementation  (ie:  16  channel 
acquisition  cards  of  sufficient  dynamic  range  and  frequency  of  operation  are  readily 
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available  now;  the  eleven  near  field  array  would  require  only  a  single  eard  for 
implementation.) 

C.  SUMMARY 

It  was  desired  to  verily  the  teehnique  with  aetual  measurement  data,  however,  the 
present  eonliguration  of  the  Aeoustie  Pressure  Tank  Faeility’s  measurement  system  (  4 
channel  ;  VXI-based)  precluded  a  reasonably  fast  or  convenient  way  to  measure  the 
complex  response  of  an  eleven  element  near  field  array.  This  latter  array  also  does  not 
presently  exist,  but  likely  could  have  been  conjured  up  within  existing  USRD  inventories 
if  so  required.  It  is  anticipated  that  this  technique  will  be  looked  at  with  an  updated 
acquisition  system  and  perhaps  a  specially  designed  array  at  some  future  date. 

An  error  study  was  undertaken  with  the  simulated  data  sets.  This  error  study 
examined  errors  in  range,  vertical  placement  and  angular  orientation  of  the  near  field 
array.  This  study  confirmed  that  the  technique  (at  least  for  the  given  set-up)  was 
extremely  robust  and  did  not  require  exacting  tolerances  to  achieve  qualifiable  results. 
This  was  especially  gratifying  for  the  test  environment  in  which  the  technique  is 
ultimately  intended  to  be  used. 

The  technique  will  not,  however,  be  appropriate  for  all  devices  and  all  operating 
frequencies.  The  reliance  on  COTS  computing  hardware  and  commercial  software 
products  does  limit  the  size  of  the  matrices  that  can  be  reasonably  handled.  As 
mentioned  previously,  the  solution  set  example  used  a  relatively  low  frequency  of  about 
4775  Hz  but  yet  still  required  the  solution  of  30,240  “unknowns”  for  its  15,120  individual 
patches. 

In  summary,  the  technique  does  present  itself  favorably  and  is  highly  likely  to  be 
looked  at  for  future  implementation  at  USRD. 
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FAR  FIELD  EXTRAPOLATION  TECHNIQUE  USING 
CHIEF  ENCLOSING  SPHERE  DEDUCED  PRESSURES  AND  VELOCITIES 

FIGURES 
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Figure  7  -  Enclosing  Sphere  Section  (numv  =84;  symmetry=180) 


NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  8  -  Theoretical  Near  Field  Response  of  Finear  Array 
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FAR  FIELD  DATA  with  k  =  20 :  range  =  11 30.9734  meters;  array  length  =  1 .885  meters 


Figure  9  -  Theoretical  Far  Field  Response  of  Linear  Array 
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Figure  10  -  Near  Field  Flydrophone  Difference  Data 
alpha0=l;  betaO  =  1 


Ma-new  -  min(HYDROPHONE),  max(ZeroO),  max(Zero1);  alpha0  =  1  ;  alpha  =  1.9789  ;  betaO  =  1  ;  beta  =  998.5472 


Figure  1 1  -  Ma-new  Components 
alpha0=l;  betaO  =  1 
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Ma-new  [  min(HYDROPHOI\IE)  -  max(ZeroO)  ]  (dB) ;  alphaO  =  1  ;  betaO  =  1 


Figure  12  -  Ma-new  [min(Hydrophone)  -  max(ZeroO)]  Data 
alphaO=l;  betaO  =  1 


Ma-new  [  min(HYDROPHONE)  -  max(Zero1)  ]  (dB) ;  alphaO  =  1  ;  betaO  =  1 


Figure  13  -  Ma-new  [min(Hydrophone)  -  max(Zerol)]  Data 
alphaO=l;  betaO  =  1 
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Figure  14  -  Ma-new  and  Ma  (original)  Data  Sets 
alphaO=l;  betaO  =  1 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  15  -  Near  Field  Data  -  computed  and  simulated 
alpha0=l;  betaO  =  1 


NF  DATA  Difference  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters 


Figure  16  -  Near  Field  Data  Difference 
alphaO=l;  betaO  =  1 


59 


FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  17  -  Far  Field  Data  -  computed  and  simulated 
alpha0=l;  betaO  =  1 


FF  DATA  Difference(on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters 
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Figure  18  -  Far  Field  Data  Difference 
alpha0=l;  betaO  =  1 
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Figure  19  -  Near  Field  Flydrophone  Difference  Data 
alphaO=l;  betaO  =10 


Ma-new-  min(HYDROPHONE),  max(ZeroO),  max(Zero1);  alphaO  =  1  ;  alpha  =  1.9789  ;  betaO  =  10  ;  beta  =  9985.4718 


Figure  20  -  Ma-new  Components 
alpha0=l;  betaO  =10 
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Ma-new  [  min(HYDROPHONE)  -  max(ZeroO)  |  (dB) ;  alphaO  =  1  ;  betaO  =  10 


Figure  21  -  Ma-new  [min(Hydrophone)  -  max(ZeroO)]  Data 
alpha0=l;  betaO  =10 


Ma-new  [  min(HYDROPHONE)  -  max(Zero1)  ]  (dB) ;  alphaO  =  1  ;  betaO  =  10 


Figure  22  -  Ma-new  [min(Hydrophone)  -  max(Zerol)]  Data 
alpha0=l;  betaO  =10 
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Figure  23  -  Ma-new  and  Ma  (original)  Data  Sets 
alphaO=l;  betaO  =10 
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HYDROPHONE  DATA  -  dB[(compufecl-actual)l ;  alphaO  =  1  ;  betaO  =  100 


Figure  24  -  Near  Field  Flydrophone  Difference  Data 
alpha0=l;  betaO  =  100 


Ma-new-  min(HYDROPHONE),  max(ZeroO),  max(Zero1);  alphaO  =  1  ;  alpha  =  1 .9789  ;  betaO  =  100  ;  beta  =  99854.7185 


Figure  25  -  Ma-new  Components 
alpha0=l;  betaO  =  100 
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Figure  26  -  Ma-new  [min(Hydrophone)  -  max(ZeroO)]  Data 
alphaO=l;  betaO  =  100 
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Figure  27  -  Ma-new  [min(Hydrophone)  -  max(Zerol)]  Data 
alphaO=l;  betaO  =  100 


400 
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Figure  28  -  Ma-new  and  Ma  (original)  Data  Sets 


alphaO=l;  betaO  =  100 


66 


HYDROPHONE  DATA  -  dB[(compufecl-actual)] ;  alphaO  =  1  ;  betaO  =  500 


Figure  29  -  Near  Field  Flydrophone  Difference  Data 
alpha0=l;  betaO  =  500 


Ma-new-  min(HYDROPHONE),  max(ZeroO),  max(Zero1):  alphaO  =  1  ;  alpha  =  1.9789  ;  betaO  =  500  ;  beta  =  499273.5925 


Figure  30  -  Ma-new  Components 
alpha0=l;  betaO  =  500 
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Figure  31  -  Ma-new  [min(Hydrophone)  -  max(ZeroO)]  Data 
alphaO=l;  betaO  =  500 


Ma-new  [  min(HYDROPHONE)  -  max(Zero1)  ]  (dB) ;  alphaO  =  1  ;  betaO  =  500 


Figure  32  -  Ma-new  [min(Hydrophone)  -  max(Zerol)]  Data 
alpha0=l;  betaO  =  500 
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Figure  33  -  Ma-new  and  Ma  (original)  Data  Sets 


alphaO=l;  betaO  =  500 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  34  -  Near  Field  Data  -  computed  and  simulated 
alpha0=l;  betaO  =  500 
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Figure  35  -  Near  Field  Data  Difference 
alphaO=l;  betaO  =  500 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  36  -  Far  Field  Data  -  computed  and  simulated 
alpha0=l;  betaO  =  500 


Figure  37  -  Far  Field  Data  Difference 
alphaO=l;  betaO  =  500 


71 


X  10°  HYDROPHONE  DATA  -  dB[(computed-actual)] ;  alphaO  =  1  ;  betaO  =  1000 


Figure  38  -  Near  Field  Flydrophone  Difference  Data 
alpha0=l;  betaO  =  1000 


Ma-new  -  from  pv  solution 


Ma  -  original 


50  100  150 


50  100  150 


Figure  39  -  Ma-new  and  Ma  (original)  Data  Sets 
alpha0=l;  betaO  =  1000 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  40  -  Near  Field  Data  -  computed  and  simulated 
alpha0=l;  betaO  =  1000 
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Figure  41  -  Near  Field  Data  Difference 
alphaO=l;  betaO  =  1000 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  42  -  Far  Field  Data  -  computed  and  simulated 
alpha0=l;  betaO  =  1000 


FF  DATA  Difference(on  axis)  with  k  =  20 ;  range  =  1130.9734  meters 


Figure  43  -  Far  Field  Data  Difference 
alphaO=l;  betaO  =  1000 
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Figure  44  -  Near  Field  Flydrophone  Difference  Data 
alpha0=10;  betaO  =  1 


Ma-new-  min(HYDROPHONE),  max(ZeroO),  max(Zero1);  alphaO  =  10  ;  alpha  =  19.7892  ;  betaO  =  1  ;  beta  =  998.5472 


Figure  45  -  Ma-new  Components 
alpha0=10;  betaO  =  1 
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Ma-new  [  min(HYDROPHONE)  -  max(ZeroO)  |  (dB) ;  alphaO  =  10  ;  betaO  =  1 


Figure  46  -  Ma-new  [min(Hydrophone)  -  max(ZeroO)]  Data 
alpha0=10;  betaO  =  1 


Ma-new  [  min(HYDROPHONE)  -  max(Zero1)  |  (dB) ;  alphaO  =  10  ;  betaO  =  1 


Figure  47  -  Ma-new  [min(Hydrophone)  -  max(Zerol)]  Data 
alpha0=10;  betaO  =  1 
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Figure  48  -  Ma-new  and  Ma  (original)  Data  Sets 


alpha0=10;  betaO  =  1 
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Figure  49  -  Near  Field  Flydrophone  Difference  Data 
alpha0=100;  betaO  =  1 


Ma-new-  min(HYDROPHONE),  max(ZeroO),  max(Zero1);  alphaO  =  100  ;  alpha  =  197.8919  ;  betaO  =  1  ;  beta  =  998.5472 


Figure  50  -  Ma-new  Components 
alpha0=100;  betaO  =  1 
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Ma-new  [  min(HYDROPHONE)  -  max(ZeroO)  ]  (dB) ;  alphaO  =  100  ;  betaO  =  1 


Figure  51  -  Ma-new  [min(Hydrophone)  -  max(ZeroO)]  Data 
alpha0=100;  betaO  =  1 


Ma-new  [  min(HYDROPHONE)  -  max(Zero1)  ]  (dB) ;  alphaO  =  100  ;  betaO  =  1 


Figure  52  -  Ma-new  [min(Hydrophone)  -  max(Zerol)]  Data 
alpha0=100;  betaO  =  1 
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Figure  53  -  Ma-new  and  Ma  (original)  Data  Sets 


alpha0=100;  betaO  =  1 
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Figure  54  -  Near  Field  Flydrophone  Difference  Data 
alpha0=500;  betaO  =  1 


Ma-new-  min(HYDROPHONE),  max(ZeroO),  max(Zero1);  alphaO  =  500  ;  alpha  =  989.4593  ;  betaO  =  1  ;  beta  =  998.5472 


Figure  55  -  Ma-new  Components 
alpha0=500;  betaO  =  1 
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Ma-new  [  min(HYDROPHONE)  -  max(ZeroO)  ]  (dB) ;  alphaO  =  500  ;  betaO  =  1 
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Figure  56  -  Ma-new  [min(Hydrophone)  -  max(ZeroO)]  Data 
alpha0=500;  betaO  =  1 
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Figure  57  -  Ma-new  [min(Hydrophone)  -  max(Zerol)]  Data 
alpha0=500;  betaO  =  1 
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Figure  58  -  Ma-new  and  Ma  (original)  Data  Sets 


alpha0=500;  betaO  =  1 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  59  -  Near  Field  Data  -  computed  and  simulated 
alpha0=500;  betaO  =  1 
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Figure  60  -  Near  Field  Data  Difference 
alpha0=500;  betaO  =  1 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  61  -  Far  Field  Data  -  computed  and  simulated 
alpha0=500;  betaO  =  1 


FF  DATA  Difference(on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters 


Figure  62  -  Far  Field  Data  Difference 
alpha0=500;  betaO  =  1 
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HYDROPHONE  DATA  -  dB[(computed-actual)] ;  alphaO=  1  ;  betaO=  1 


Figure  63  -  Near  Field  Flydrophone  Difference  Data 
alphaO=l;  betaO=l;  range  error  =  -2  *  lambda 


Figure  64  -  Ma-new  and  Ma  (original)  Data  Sets 
alphaO^l;  betaO^l;  range  error  =  -2  *  lambda 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  65  -  Near  Field  Data  -  eomputed  and  simulated 
alpha0=l;  betaO^l;  range  error  =  -2  *  lambda 


X 10'^  NF  DATA  Difference  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters 
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Figure  66  -  Near  Field  Data  Differenee 
alpha0=l;  beta0=l;  range  error  =  -2  *  lambda 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  67  -  Far  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  range  error  =  -2  *  lambda 


FF  DATA  Difference(on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters 


Figure  68  -  Far  Field  Data  Difference 
alphaO^l;  betaO^l;  range  error  =  -2  *  lambda 
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HYDROPHONE  DATA  -  dB[(computed-actual)] ;  alphaO=  1  ;  betaO=  1 


Figure  69  -  Near  Field  Flydrophone  Difference  Data 
alphaO=l;  betaO=l;  range  error  =  +2  *  lambda 


Figure  70  -  Ma-new  and  Ma  (original)  Data  Sets 
alpha0=l;  beta0=l;  range  error  =  +2  *  lambda 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  71  -  Near  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  range  error  =  +2  *  lambda 


NF  DATA  Difference  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters 


Figure  72  -  Near  Field  Data  Difference 
alpha0=l;  beta0=l;  range  error  =  +2  *  lambda 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  73  -  Far  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  range  error  =  +2  *  lambda 
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Figure  74  -  Far  Field  Data  Difference 
alpha0=l;  beta0=l;  range  error  =  +2  *  lambda 
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HYDROPHONE  DATA  -  dB[(computed-actual)] ;  alphaO=  1  ;  betaO=  1 


Figure  75  -  Near  Field  Flydrophone  Difference  Data 
alphaO=l;  betaO=l;  range  error  =  -1  cm 


Figure  76  -  Ma-new  and  Ma  (original)  Data  Sets 
alphaO=l;  betaO=l;  range  error  =  -1  cm 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  77  -  Near  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  range  error  =  -1  cm 


X 10'^  NF  DATA  Difference  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters 
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Figure  78  -  Near  Field  Data  Difference 
alpha0=l;  beta0=l;  range  error  =  -1  cm 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  79  -  Far  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  range  error  =  -1  cm 


FF  DATA  Difference(on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters 


Figure  80  -  Far  Field  Data  Difference 
alpha0=l;  beta0=l;  range  error  =  -1  cm 
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Figure  81  -  Near  Field  Flydrophone  Difference  Data 
alpha0=l;  beta0=l;  range  error  =  +1  cm 


Figure  82  -  Ma-new  and  Ma  (original)  Data  Sets 
alpha0=l;  beta0=l;  range  error  =  +1  cm 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  83  -  Near  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  range  error  =  +1  cm 


NF  DATA  Difference  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters 


Figure  84  -  Near  Field  Data  Difference 
alpha0=l;  beta0=l;  range  error  =  +1  cm 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  85  -  Far  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  range  error  =  +1  cm 


FF  DATA  Difference(on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters 


Figure  86  -  Far  Field  Data  Difference 
alpha0=l;  beta0=l;  range  error  =  +1  cm 
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HYDROPHONE  DATA  -  dB[(computed-actual)] ;  alphaO=  1  ;  betaO=  1 


Figure  87  -  Near  Field  Flydrophone  Difference  Data 
alphaO=l;  betaO=l;  vertical  receive  array  error  =  +lambda/2 


Figure  88  -  Ma-new  and  Ma  (original)  Data  Sets 


alphaO=l;  betaO=l;  vertical  receive  array  error  =  +lambda/2 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  89  -  Near  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  +lambda/2 


X 10'^  NF  DATA  Difference  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters 
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Figure  90  -  Near  Field  Data  Difference 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  +lambda/2 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  91  -  Far  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  +lambda/2 
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Figure  92  -  Far  Field  Data  Difference 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  +lambda/2 
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Figure  93  -  Near  Field  Flydrophone  Difference  Data 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  +lambda 


Figure  94  -  Ma-new  and  Ma  (original)  Data  Sets 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  +lambda 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  95  -  Near  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  +lambda 


NF  DATA  Difference  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters 


Figure  96  -  Near  Field  Data  Difference 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  +lambda 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  97  -  Far  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  +lambda 


FF  DATA  Difference(on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters 


Figure  98  -  Far  Field  Data  Difference 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  +lambda 
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Figure  99  -  Near  Field  Flydrophone  Difference  Data 
alphaO=l;  betaO=l;  vertical  receive  array  error  =  -lambda/2 


Figure  100  -  Ma-newand  Ma  (original)  Data  Sets 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  -lambda/2 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  101  -  Near  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  -lambda/2 
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Figure  102  -  Near  Field  Data  Difference 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  -lambda/2 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 


Figure  103  -  Far  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  -lambda/2 
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Figure  104  -  Far  Field  Data  Difference 
alpha0=l;  beta0=l;  vertical  receive  array  error  =  -lambda/2 
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Figure  105  -  Near  Field  Flydrophone  Difference  Data 
alpha0=l;  beta0=l;  azimuthal  error  =  +10  degrees 


Figure  106  -  Ma-newand  Ma  (original)  Data  Sets 
alpha0=l;  beta0=l;  azimuthal  error  =  +10  degrees 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  5.3541  meters;  array  length  =  1 .885  meters 


Figure  107  -  Near  Field  Data  -  computed  and  simulated 
alpha0=l;  beta0=l;  azimuthal  error  =  +10  degrees 


Figure  108  -  Near  Field  Data  Difference 
alpha0=l;  beta0=l;  azimuthal  error  =  +10  degrees 
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FAR  FIELD  DATA  (on  axis)  with  k  =  20 :  range  =  1130.9734  meters;  array  length  =  1 .885  meters 
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Figure  109  -  Far  Field  Data  -  computed  and  simulated 
alphaO=l;  betaO=l;  azimuthal  error  =  +10  degrees 


FF  DATA  Difference(on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters 


Figure  1 10  -  Far  Field  Data  Difference 
alpha0=l;  beta0=l;  azimuthal  error  =  +10  degrees 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  2.0817  meters;  array  length  =  1 .885  meters 


Figure  1 1 1  -  Near  Field  Data  -  computed  and  simulated 
Using  2.082  meter  range,  1.571  meter  long,  1 1  element  near  field  array  data 


FAR  FIELD  DATA  (on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters;  array  length  =  1 .885  meters 


Figure  1 12  -  Far  Field  Data  -  computed  and  simulated 
Using  2.082  meter  range,  1.571  meter  long,  1 1  element  near  field  array  data 


no 


Figure  1 13  -  Ma-new  and  Ma  (original)  Data  Sets 
Using  2.082  meter  range,  1.571  meter  long,  1 1  element  near  field  array  data 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  1 4658  meters;  array  length  =  1 .885  meters 


Figure  1 14  -  Near  Field  Data  -  computed  and  simulated 
Using  1.466  meter  range,  1.571  meter  long,  1 1  element  near  field  array  data 


FAR  FIELD  DATA  (on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters;  array  length  =  1 .885  meters 


Figure  1 15  -  Far  Field  Data  -  computed  and  simulated 
Using  1.466  meter  range,  1.571  meter  long,  1 1  element  near  field  array  data 
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Figure  1 16  -  Ma-new  and  Ma  (original)  Data  Sets 
Using  1.466  meter  range,  1.571  meter  long,  1 1  element  near  field  array  data 
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NEAR  FIELD  DATA  (Rec  Ele  #6)  with  k  =  20  ;  range  =  1 4658  meters;  array  length  =  1 .885  meters 


Figure  1 17  -  Near  Field  Data  -  computed  and  simulated 
Using  1.466  meter  range,  1.571  meter  long,  1 1  element  near  field  array  data 


FAR  FIELD  DATA  (on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters;  array  length  =  1 .885  meters 


Figure  1 18  -  Far  Field  Data  -  computed  and  simulated 
Using  1.466  meter  range,  1.571  meter  long,  1 1  element  near  field  array  data 
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Figure  1 19  -  Ma-new  and  Ma  (original)  Data  Sets 


Using  1.466  meter  range,  1.571  meter  long,  1 1  element  near  field  array  data 
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NEAR  FIELD  DATA  (Rec  Ele  #8)  with  k  =  20  ;  range  =  1 4658  meters;  array  length  =  1 .885  meters 


Figure  120  -  Near  Field  Data  -  computed  and  simulated 
Using  1.466  meter  range,  2.199  meter  long,  15  element  near  field  array  data 


FAR  FIELD  DATA  (on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters;  array  length  =  1 .885  meters 


Figure  121  -  Far  Field  Data  -  computed  and  simulated 
Using  1.466  meter  range,  2.199  meter  long,  15  element  near  field  array  data 
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Figure  122  -  Ma-newand  Ma  (original)  Data  Sets 


Using  1.466  meter  range,  2.199  meter  long,  15  element  near  field  array  data 
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NEAR  FIELD  DATA  (Rec  Ele  #7)  with  k  =  20  :  range  =  1 .885  meters;  array  length  =  1 .885  meters 


Figure  123  -  Near  Field  Data  -  computed  and  simulated 
Using  1.885  meter  range,  1.885  meter  long,  13  element  near  field  array  data 


FAR  FIELD  DATA  (on  axis)  with  k  =  20 ;  range  =  1 130.9734  meters;  array  length  =  1 .885  meters 


Figure  124  -  Far  Field  Data  -  computed  and  simulated 
Using  1.885  meter  range,  1.885  meter  long,  13  element  near  field  array  data 
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Figure  125  -  Ma-newand  Ma  (original)  Data  Sets 
Using  1.885  meter  range,  1.885  meter  long,  13  element  near  field  array  data 
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APPENDIX  A  -  EXAMPLE  AB  MATRIX 


The  following  pages  show  seetions  of  a  full  augmented  AB_Matrix  as  identified 
in  this  effort.  This  example  is  provided  to  assist  in  visualizing  the  matrix  eomponents, 
their  relationship  to  one  another  and  the  overall  eirculant  nature  of  the  matrix.  The 
example  provided  details  an  enclosing  sphere  of  very  small  order.  Due  to  this  small  size, 
this  matrix  would  not  be  suited  for  actual  use  but  makes  it  ideal  for  instructional 
purposes. 

For  this  example,  an  enclosing  sphere  was  constructed  with  the  following 
dimensions; 

numv  =  5  this  represents  the  number  of  patches  in  the  “north-south”  direction 

sph_sym  =10  this  represents  the  number  of  “slices”  that  comprise  the  sphere 

numh  =  3  this  represents  the  number  of  elements  in  the  near  field  linear  array 

The  enclosing  sphere  so  obtained  would  look  like  the  following; 


121 


As  discussed  in  the  write-up,  the  eonstituents  of  the  AB_Matrix  are  of  the 
defined  sizes  listed  below: 

Matrix 

A 
B 

delG 
G 

PE 
VE 

where: 

PE  =  Pressure  equivalence  matrix 
VE  =  Veloeity  equivalenee  matrix 

The  AB_Matrix  is  formed  from  the  eoneatenation  of  the  above  eomponents. 
This  matrix  takes  the  form  of: 


Ai 

A2 

B2 

Aso 

Bso 

delGi 

Gi 

delG2 

G2 

delGso  Gso 

PEi 

0 

PE2 

0 

PEso 

0 

0 

VEi 

0 

VE2 

0 

VE50 

The  size  of  one  “slice”  of  the  AB_matrix  is  represented  by  the  above 
arrangement.  Eor  the  given  example  parameters,  this  matrix  is  of  size  16  x  100.  Since 
the  full  matrix  is  cireulant  due  to  the  symmetry  of  the  ehosen  enelosing  surfaee  (a 
sphere),  the  full  AB_matrix  is  of  size  160  x  100.  (A  matrix  is  eonsidered  to  be  eireulant 


Size  (rows  x  columns) 

Example  Matrix 

numv  X  (sph  numv  *  sph  sym) 

(4  it. 

5  x50 

5x50 

numh  X  (sph  numv  *  sph  sym) 

3  X  50 

3  X  50 

(  sph  numv-l)  x  (sph  numv  *  sph  sym) 

a  a 

4x50 

4x50 
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if  each  of  its  rows,  or  columns,  is  a  circular  shift  of  the  previous  row,  or  column.  If  the 
matrix  elements  are  considered  to  be  matrices  themselves  and  they  obey  the  same 
circulant  behavior  as  identified  above,  the  matrix  is  more  correctly  identified  as  being 
block  circulant.) 

The  pages  that  follow  show  some  sections  of  the  full  example  AB_matrix.  It  was 
not  practical  nor  of  any  real  value  to  include  the  entire  16,000  element  matrix.  Rather, 
some  representative  subset  was  provided.  The  matrix  has  been  numbered  for  reference 
purposes  with  column  numbers  appearing  across  the  top  or  bottom  and  row  numbers 
appearing  on  the  left  hand  side. 

Since  the  array  is  circulant,  the  elements  defined  in  the  first  slice  (Rows  1  thru  16) 
and  (Columns  1  thru  100)  will  repeat  in  a  predictable  manner  throughout  the  rest  of  the 
matrix.  Rows  1  thru  16  in  Column  1  have  been  highlighted  in  red  to  show  this  repetition 
in  subsequent  slices.  The  data  that  appears  in  (Rows  1  thru  16),  column  1  also  appears  in 
(Rows  17  thru  32),  columnll  and  will  appear  in  (Rows  33  thru  48),  Column  21  and  in 
other  slices  up  until  its  final  appearance  in  (Rows  145  thru  160),  Column  91. 

Examination  of  the  other  matrix  elements  shows  that  subsequent  slices  are  created 
by  shifting  the  present  slice  sph  sym  (in  our  case  this  is  10)  positions  to  the  right  and 
allowing  those  columns  shifted  past  the  last  column  (in  our  case  this  is  100)  to  be 
wrapped  back  to  the  beginning  of  the  slice.  This  occurs  a  total  of  (sph_sym-l)  times 
yielding  the  final  matrix. 

Another  significant  point  in  the  construction  of  this  matrix  is  the  manner  in  which 
the  adjacent  patch  pressure  equivalence  (PE)  and  velocity  equivalence  (VE)  is 
represented.  Along  any  single  row  of  PE  or  VE,  there  will  appear  a  single  1  value 
followed  immediately  by  a  single  -1  value;  all  other  elements  in  these  rows  are  zero. 
This  forms  the  conditions  of  (pi  -  pi+i)  or  (vi  -  Vi+i)  that  is  required  for  adjacent  patch 
locations  on  the  enclosing  sphere.  An  examination  of  the  matrix  will  actually  show  1,0,- 
1  (as  in  row  9,  columns  1  thru  3)  that  is  used  to  implement  this  condition.  The  zero  in 
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this  case  is  actually  associated  with  the  velocity  term  and  thus  does  not  interrupt  the 
adjacency  requirement  between  successive  pressure  terms.  This  is  an  artifact  of  the 
interleaving  of  the  A  and  B  matrices. 

The  AB_matrix  represents  the  known  terms  in  the  overdetermined  system  that  it 
creates.  That  is,  in  the  example  provided,  there  are  100  total  unknowns.  These 
correspond  to  50  unknown  pressures  and  50  unknown  velocities.  A  single  to-be-solved- 
for  pressure  that  will  multiply  all  elements  of  column  1  will  be  generated  via  the  least 
squares  solution  technique.  A  single  to-be-solved-  for  velocity  that  will  multiply  all 
elements  of  column  2  will  also  be  generated.  Thus  pressure  values  corresponding  to  all 
odd  numbered  columns  and  velocity  values  corresponding  to  all  even  numbered  columns 
are  ultimately  generated. 

Before  leaving  this  appendix,  it  is  significant  to  provide  some  additional 
comments  on  the  symmetry  requirement  of  the  aforementioned  components.  In  the 
MATLAB  implementation,  a  single  slice  of  the  augmented  AB_matrix  is  formed.  Its 
symmetry  conditions  allow  either  the  full  AB_matrix  to  be  determined  (as  was  illustrated 
in  this  appendix)  or  it  may  be  used  directly  with  “rotational  solve”  tools  available  in 
CHIEF  for  MATLAB  to  obtain  a  solution.  This  latter  usage  is  preferred  since  it  is 
computationally  more  efficient  and  imposes  significantly  less  memory  demands  than  the 
full  defined  AB_matrix. 

To  utilize  these  “rotational  solve”  tools,  however,  requires  that  the  AB_matrix 
slice  be  representative  of  a  section  of  the  full  block  circulant  AB_matrix.  It  is  for  this 
latter  reason  that  such  care  was  given  to  the  definition  of  the  enclosing  sphere,  the 
placement  of  the  near  field  element  points  and  the  manner  in  which  the  adjacency 
conditions  were  defined.  The  symmetrical  enclosing  sphere  satisfies  the  basic  required 
circulant  behavior  since  it  is  related  to  the  field  points  in  a  very  predictable  manner.  That 
is,  a  rotation  of  the  enclosing  sphere  by  one  slice  (with  sph  sym  being  the  total  number 
of  increments  in  a  complete  rotation)  yields  the  same  values  of  A,  B  ,  G  and  delG  with 
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the  exception  that  their  placement  in  the  matrix  has  shifted.  In  the  example  of  this 
appendix,  the  columns  of  data  shift  10  columns  to  the  right  for  each  incremental  rotation. 
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Augmented 
AB  Matrix 


1 

2 

3 

4 

5 

6 

A  and  deIG 

B  and  G 

A  and  deIG 

B  and  G 

A  and  deIG 

B  and  G 

1 

-2.995  +  0.173761 

3.5684  +  0.921 34i 

-0.093182  +  0.151741 

0.54411  -0.19015i 

-0.013249-0.358891 

-0.032079  +  0.62409i 

2 

-0.10782  +  0.00369071 

0.52906  +  0.11537i 

-2.6184  +  0.14508i 

1.0352  -0.82771  i 

-0.11272  +  0.20776i 

0.7104-0.440971 

3 

0.017724-0.102391 

-0.047656  +  0.21 304i 

-0.12407  +  0.16402i 

0.69458- 0.36421  i 

-2.6275  -  0.083952i 

0.67693  +  0.59722i 

4 

0.023804  +  0.097903i 

-0.01765 -0.1466i 

-0.02217 -0.28557i 

0.00057447  +  0.51 634i 

-0.11272  +  0.20776i 

0.7104 -0.44097i 

5 

0.17511  +0.27727i 

-0.17103 -0.30889i 

-0.10198  +  0.0442711 

0.11917 -0.098908i 

-0.013249 -0.35889i 

-0.032079  +  0.62409i 

6 

4.6858  + 2.4081  i 

-10.891  -  14.991i 

1.2014 -6.6544i 

-4.6272  +  28.753i 

2.9742  -  6.3597i 

-26.025  +  12.142i 

7 

-3.0872  +  1.95491 

17.786  +  3.1931i 

4.5592  +  0.99307i 

-29.777  -  6.7274i 

6.4483- 3.51 3i 

-29.578  -  5.4269i 

8 

-2.5756-0.0910331 

3.9504  +  17.378i 

1.8408+  1.0087i 

-23.936-  17.67i 

2.9742  -  6.3597i 

-26.025  +  12.142i 

9 

1 

0 

-1 

0 

0 

0 

10 

0 

0 

1 

0 

-1 

0 

11 

0 

0 

0 

0 

1 

0 

12 

0 

0 

0 

0 

0 

0 

13 

0 

1 

0 

-1 

0 

0 

14 

0 

0 

0 

1 

0 

-1 

15 

0 

0 

0 

0 

0 

1 

16 

0 

0 

0 

0 

0 

0 

17 

0.12274 -0.024667i 

-0.91857-0.719031 

0.12452  +  0.129671 

-0.15791  -0.45878i 

-0.25971  -0.19699i 

0.40766  +  0.35186i 

18 

-0.0027186-0. 048651  i 

-0.039756  + 0.0831 06i 

-0.1168-0.163941 

0.45809  +  0.83399i 

0.05291  -0.014402i 

-0.1824-0.0310351 

19 

-0.028833  +  0.041 861  i 

0.079527  -  0.06373i 

0.13329 -0.062714i 

-0.34442  +  0.075509i 

-0.11712  +  0.22948i 

0.70567- 0.51 81 5i 

20 

0.14422 -0.0026827i 

-0.18576-0.0217151 

-0.12969  +  0.00516171 

0.20041  -  0.0069979i 

0.05291  -0.014402i 

-0.1824-0.0310351 

21 

0.25895  + 0.24221  i 

-0.26248  -  0.27254i 

-0.033291  +  0.026667i 

0.0081437  -  0.066227i 

-0.25971  -0.19699i 

0.40766  +  0.35186i 

22 

-3.4698  +  0.78644i 

8.1228  +  3.2625i 

-7.4585  + 0.94841  i 

10.618 -6.5525i 

-18.52  +  17.577i 

24.984- 26.491  i 

23 

0.38449  -  2.2907i 

-5.7626  +  2.9136i 

-3.2486- 5.71 39i 

9.1501  +9.1485i 

-26.567  +  4.8747i 

37.351  -  9.0999i 

24 

1.2807-  1.1 678i 

-3.6595  -  2.6948i 

-2.6672-5.02121 

6.7898+  10.278i 

-18.52  +  17.577i 

24.984- 26.491  i 

25 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

27 

0 

0 

0 

0 

0 

0 

28 

0 

0 

0 

0 

0 

0 

29 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 

31 

0 

0 

0 

0 

0 

0 

32 

0 

0 

0 

0 

0 

0 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


7 

A  and  deIG 


8 

B  and  G 


9 

A  and  deIG 


10 

B  and  G 


11 

A  and  deIG 


12 

B  and  G 


-0.10198  +  0.044271  i 

0.11917-0.0989081 

0.17511  +  0.27727i 

-0.02217-0.285571 

0.00057447  +  0.51 634i 

0.023804  +  0.097903i 

-0.12407  +  0.164021 

0.69458-0.364211 

0.017724-0.102391 

-2.6184  +  0.145081 

1.0352-0.827711 

-0.10782  +  0.00369071 

-0.093182  +  0.151741 

0.54411  -0.19015i 

-2.995  +  0.173761 

1.8408  +  1.0087i 

-23.936 -  17.67i 

-2.5756-0.0910331 

4.5592  +  0.99307i 

-29.777  -  6.7274i 

-3.0872  +  1.9549i 

1.2014-6.65441 

-4.6272  +  28.753i 

4.6858  +  2.40811 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

1 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

1 

0 

-0.033291  +  0.026667i 

0.0081437  -  0. 0662271 

0.25895  +  0.242211 

-0.12969  +  0.00516171 

0.20041  -  0.0069979i 

0.14422-0.00268271 

0.13329-0.0627141 

-0.34442  +  0.075509i 

-0.028833  +  0.041 861  i 

-0.1168-0.163941 

0.45809  +  0.83399i 

-0.0027186-0.0486511 

0.12452  +  0.129671 

-0.15791  -0.45878i 

0.12274-  0.0246671 

-2.6672-5.02121 

6.7898  +  10.278i 

1.2807-  1.16781 

-3.2486-5.71391 

9.1501  +  9.1485i 

0.38449  -  2.2907i 

-7.4585  +  0.948411 

10.618-6.55251 

-3.4698  +  0.78644i 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.17103-0.308891 

0.12274-0.0246671 

-0.91857-0.719031 

-0.01765-0.14661 

-0.0027186-0.0486511 

-0.039756  +  0.0831061 

-0.047656  +  0.213041 

-0.028833  +  0.041 861  i 

0.079527  -  0.06373i 

0.52906  +  0.115371 

0.14422-0.00268271 

-0.18576-0.0217151 

3.5684  +  0.921341 

0.25895  +  0.242211 

-0.26248  -  0.27254i 

3.9504  +  17.378i 

-1.4722-  1.147i 

0.75928  +  3.8456i 

17.786  +  3.19311 

1.7204-0.407311 

-4.0252  -  2.8299i 

-10.891  -  14.991i 

1.3635  +  0.954431 

1.3403-6.03111 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

-0.26248  -  0.27254i 

-2.995  +  0.173761 

3.5684  +  0.921341 

-0.18576-0.0217151 

-0.10782  +  0.00369071 

0.52906  +  0.115371 

0.079527  -  0.06373i 

0.017724-0.102391 

-0.047656  +  0.213041 

-0.039756  +  0.0831061 

0.023804  +  0.0979031 

-0.01765-0.14661 

-0.91857-0.719031 

0.17511  +0.277271 

-0.17103-0.308891 

-3.6595  -  2.6948i 

4.6858  +  2.40811 

-10.891  -  14.991i 

-5.7626  +  2.91361 

-3.0872  +  1.95491 

17.786  +  3.19311 

8.1228  +  3.26251 

-2.5756-0.0910331 

3.9504+  17.378i 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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19 

A  and  deIG 


20 

B  and  G 


21 

A  and  deIG 


22 

B  and  G 


23 

A  and  deIG 


24 

B  and  G 


1 

0.25895  +  0.242211 

-0.26248  -  0.27254i 

-0.016032-0.0141341 

0.044255  +  0.34347i 

-0.081084-0.159981 

0.094108  +  0.397591 

2 

0.14422-0.00268271 

-0.18576-0.0217151 

0.10769  +  0.179371 

-0.16384-0.447221 

-0.054417  +  0.0983391 

0.14141  -0.237861 

3 

-0.028833  +  0.0418611 

0.079527  -  0.06373i 

-0.013898-0.248111 

-0.016261  +0.36104i 

-0.074295-0.156541 

0.13711  +  0.224531 

4 

-0.0027186-0.0486511 

-0.039756  +  0.0831061 

-0.17666-0.267451 

0.18622  +  0.325061 

0.16616  +  0.0882261 

-0.18519-0.116331 

5 

0.12274-0.0246671 

-0.91857-0.719031 

0.41548  +  0.0498141 

-0.43229  -  0.070006i 

-0.12447-0.0801011 

0.10459  +  0.124611 

6 

1.3635  +  0.954431 

1.3403-6.03111 

0.77887  +  1.2599i 

10.297  +  4.61521 

-16.185-2.33461 

-22.008-3.81571 

7 

1.7204-0.407311 

-4.0252  -  2.8299i 

-2.9368  +  0.28087i 

-11.836  +  3.8171 

1.0893-  16.6911 

3.0386-  18.0041 

8 

-1.4722-  1.147i 

0.75928  +  3.8456i 

-3.1626-3.73351 

-9.4768  -  9.6625i 

-1.4769-  17.6951 

-2.173-  14.6661 

9 

0 

0 

0 

0 

0 

0 

10 

0 

0 

0 

0 

0 

0 

11 

0 

0 

0 

0 

0 

0 

12 

0 

0 

0 

0 

0 

0 

13 

0 

0 

0 

0 

0 

0 

14 

0 

0 

0 

0 

0 

0 

15 

0 

0 

0 

0 

0 

0 

16 

0 

0 

0 

0 

0 

0 

17 

0.17511  +0.27727i 

-0.17103-0.308891 

0.12274-0.0246671 

-0.91857-0.719031 

0.12452  +  0.129671 

-0.15791  -0.458781 

18 

0.023804  +  0.097903i 

-0.01765-0.14661 

-0.0027186-0.0486511 

-0.039756  +  0.0831 06i 

-0.1168-0.163941 

0.45809  +  0.833991 

19 

0.017724-0.102391 

-0.047656  +  0.213041 

-0.028833  +  0.041 861  i 

0.079527  -  0.06373i 

0.13329-0.0627141 

-0.34442  +  0.0755091 

20 

-0.10782  +  0.00369071 

0.52906  +  0.11537i 

0.14422-0.00268271 

-0.18576-0.0217151 

-0.12969  +  0.00516171 

0.20041  -  0.00699791 

21 

-2.995  +  0.173761 

3.5684  +  0.92134i 

0.25895  +  0.242211 

-0.26248  -  0.27254i 

-0.033291  +  0.0266671 

0.0081437  -  0.0662271 

22 

-2.5756-0.0910331 

3.9504+  17.378i 

-1.4722-  1.147i 

0.75928  +  3.8456i 

2.8948-3.31811 

3.8416-  13.271 

23 

-3.0872  +  1.9549i 

17.786  +  3.1931i 

1.7204-0.407311 

-4.0252  -  2.8299i 

5.1313  +  2.27991 

16.787  +  3.83071 

24 

4.6858  +  2.40811 

-10.891  -  14.991i 

1.3635  +  0.954431 

1.3403-6.03111 

7.4033  +  3.30461 

19.042  +  9.30511 

25 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

27 

0 

0 

0 

0 

0 

0 

28 

-1 

0 

0 

0 

0 

0 

29 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 

31 

0 

0 

0 

0 

0 

0 

32 

0 

-1 

0 

0 

0 

0 

128 


127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 


0 

0 

0.25895  + 0.24221  i 
0.14422 -0.0026827i 
-0.028833 +  0.041 861  i 
-0.0027186- 0.048651  i 
0.12274 -0.024667i 
1.2807-  1.1 678i 
0.38449  -  2.2907i 
-3.4698  +  0.78644i 
0 
0 
0 
0 
0 
0 
0 
0 

0.41548  +  0.049814i 
-0.17666 -0.26745i 
-0.013898- 0.2481 1i 
0.10769 +  0.17937i 
-0.016032 -0.014134i 
-2.8241  +  4.751 3i 
2.6929  +  6.5929i 
5.4551  -  6.5507i 
0 
0 
0 
0 
0 
0 
0 
0 
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0 

-1 

-0.26248  -  0.27254i 
-0.18576 -0.021715i 
0.079527  -  0.06373i 
-0.039756  +  0.0831061 
-0.91857-0.719031 
-3.6595  -  2.6948i 
-5.7626 +  2.9136i 
8.1228  +  3.26251 
0 
0 
0 
0 
0 
0 
0 
0 

-0.43229  -  0.070006i 
0.18622  +  0.32506i 
-0.016261  +  0.36104i 
-0.16384-0.447221 
0.044255  +  0.34347i 
6.5666-  10.689i 
-5.0602-  11. 867i 
-8.0569  +  10.1371 
0 
0 
0 
0 
0 
0 
0 
0 
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0 

0 

-2.995  + 0.1 7376i 

-0.10782  +  0.0036907! 

0.017724-0.10239! 

0.023804  +  0.097903! 

0.17511  +0.27727! 

4.6858  +  2.4081! 

-3.0872  +  1.9549! 

-2.5756  -0.091033! 

1 

0 

0 

0 

0 

0 

0 

0 

0.12274-0.024667! 

-0.0027186-0.048651! 

-0.028833  +  0.041861! 

0.14422-0.0026827! 

0.25895  +  0.24221! 

-3.4698  +  0.78644! 

0.38449  -  2.2907! 

1.2807-  1.1678! 

0 

0 

0 

0 

0 

0 

0 

0 

81 
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0 

0 

3.5684  +  0.92134! 
0.52906  +  0.11537! 
-0.047656  +  0.21304! 
-0.01765-0.1466! 
-0.17103-0.30889! 
-10.891  -  14.991! 
17.786  +  3.1931! 
3.9504+  17.378! 

0 

0 

0 

0 

1 

0 

0 

0 

-0.91857-0.71903! 
-0.039756  +  0.083106! 
0.079527  -  0.06373! 
-0.18576-0.021715! 
-0.26248  -  0.27254! 
8.1228  +  3.2625! 
-5.7626  +  2.9136! 
-3.6595  -  2.6948! 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

-0.093182  +  0.15174! 
-2.6184  +  0.14508! 
-0.12407  +  0.16402! 
-0.02217-0.28557! 
-0.10198  +  0.044271! 
1.2014  -6.6544! 
4.5592  +  0.99307! 
1.8408  +  1.0087! 

-1 

1 

0 

0 

0 

0 

0 

0 

0.12452  +  0.12967! 
-0.1168-0.16394! 
0.13329-0.062714! 
-0.12969  +  0.0051617! 
-0.033291  +  0.026667! 
-7.4585  +  0.94841! 
-3.2486-5.7139! 
-2.6672-5.0212! 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0.54411  -0.19015! 
1.0352-0.82771! 
0.69458-0.36421! 
0.00057447  +  0.51634! 
0.11917-0.098908! 
-4.6272  +  28.753! 
-29.777  -  6.7274! 
-23.936-  17.67! 

0 

0 

0 

0 

-1 

1 

0 

0 

-0.15791  -0.45878! 
0.45809  +  0.83399! 
-0.34442  +  0.075509! 
0.20041  -  0.0069979! 
0.0081437-0.066227! 
10.618-6.5525! 
9.1501  +9.1485! 
6.7898  +  10.278! 

0 

0 

0 

0 

0 

0 

0 

0 
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129 


127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 


0 

0 

0.12274 -0.024667i 
-0.0027186- 0.048651  i 
-0.028833  +  0.041 861  i 
0.14422 -0.0026827i 
0.25895  + 0.24221  i 
-1.4722-  1.147i 
1.7204- 0.40731  i 
1.3635  +  0.95443i 
0 
0 
0 
0 
0 
0 
0 
0 

-2.995  +  0.173761 
-0.10782  +  0.00369071 
0.017724-0.102391 
0.023804  +  0.097903i 
0.17511  +0.27727i 
4.6858  + 2.4081  i 
-3.0872  +  1.95491 
-2.5756-0.0910331 
1 
0 
0 
0 
0 
0 
0 
0 
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0 

0 

-0.91857-0.719031 
-0.039756  + 0.0831 06i 
0.079527  -  0.06373i 
-0.18576 -0.021715i 
-0.26248  -  0.27254i 
0.75928  +  3.8456i 
-4.0252  -  2.8299i 
1.3403- 6.031 1i 
0 
0 
0 
0 
0 
0 
0 
0 

3.5684 +  0.92134i 
0.52906  +  0.11537i 
-0.047656  +  0.21 304i 
-0.01765 -0.1466i 
-0.17103 -0.30889i 
-10.891  -  14.991i 
17.786  +  3.1931i 
3.9504+  17.378i 
0 
0 
0 
0 
1 
0 
0 
0 
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B  and  G 


0 

0 

0.12452  +  0.129671 
-0.1168-0.163941 
0.13329-0.0627141 
-0.12969  +  0.00516171 
-0.033291  +  0. 0266671 
2.8948 -3.3181i 
5.1313  +  2.27991 
7.4033  +  3.3046i 
0 
0 
0 
0 
0 
0 
0 
0 

-0.093182  +  0.151741 
-2.6184  +  0.14508i 
-0.12407  +  0.16402i 
-0.02217 -0.28557i 
-0.10198  +  0.044271  i 
1.2014 -6.6544i 
4.5592  +  0.99307i 
1.8408+  1.0087i 
-1 
1 
0 
0 
0 
0 
0 
0 
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0 

0 

-0.15791  -0.45878i 
0.45809  +  0.83399i 
-0.34442  +  0.075509i 
0.20041  -  0.0069979i 
0.0081437  -  0.066227i 
3.8416-  13.27i 
16.787 +  3.8307i 
19.042  +  9.3051i 
0 
0 
0 
0 
0 
0 
0 
0 

0.54411  -0.19015i 
1.0352  -0.82771  i 
0.69458- 0.36421  i 
0.00057447  +  0.516341 
0.11917 -0.098908i 
-4.6272  +  28.753i 
-29.777  -  6.7274i 
-23.936-  17.67i 
0 
0 
0 
0 
-1 
1 
0 
0 

94 
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0 

0 

-0.25971  -0.19699i 
0.05291  -  0.014402i 
-0.11712  +  0.22948i 
0.05291  -  0.014402i 
-0.25971  -0.19699i 
7.2473  + 21. 658i 
-8.2609  +  23.1111 
7.2473  +  21.6581 
0 
0 
0 
0 
0 
0 
0 
0 

-0.013249 -0.35889i 
-0.11272  +  0.20776i 
-2.6275  -  0. 0839521 
-0.11272  +  0.20776i 
-0.013249 -0.35889i 
2.9742  -  6.3597i 
6.4483- 3.51 3i 
2.9742  -  6.3597i 
0 
-1 
1 
0 
0 
0 
0 
0 

95 
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0 

0 

0.40766  +  0.35186i 
-0.1824-0.0310351 
0.70567-0.518151 
-0.1824-0.0310351 
0.40766  +  0.35186i 
18.329  +  37.397i 
-9.7121  +43.426i 
18.329  +  37.3971 
0 
0 
0 
0 
0 
0 
0 
0 

-0.032079  +  0.62409i 
0.7104-0.440971 
0.67693  +  0.59722i 
0.7104-0.440971 
-0.032079  +  0.62409i 
-26.025  +  12.142i 
-29.578  -  5.4269i 
-26.025  +  12.142i 
0 
0 
0 
0 
0 
-1 
1 
0 
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130 


127 

0 

0 

0 

0 

128 

0 

0 

0 

0 

129 

-0.033291  +  0.026667i 

0.0081437  -  0.066227i 

0.25895  +  0.242211 

-0.26248  -  0.27254i 

130 

-0.12969  +  0.0051617i 

0.20041  -  0.0069979i 

0.14422 -0.0026827i 

-0.18576 -0.021 71 5i 

131 

0.13329 -0.062714i 

-0.34442  +  0.075509i 

-0.028833  +  0.041 861  i 

0.079527  -  0.06373i 

132 

-0.1168- 0.1 6394i 

0.45809  +  0.83399i 

-0.0027186-0.0486511 

-0.039756  +  0.0831061 

133 

0.12452  +  0.12967i 

-0.15791  -0.45878i 

0.12274-0.0246671 

-0.91857 -0.71903i 

134 

7.4033  +  3.3046i 

19.042  +  9.3051  i 

1.3635  +  0.954431 

1.3403- 6.031 1i 

135 

5.1313  +  2.27991 

16.787  +  3.8307i 

1.7204- 0.40731  i 

-4.0252  -  2.8299i 

136 

2.8948-3.31811 

3.8416-  13.27i 

-1.4722-  1.147i 

0.75928  +  3.8456i 

137 

0 

0 

0 

0 

138 

0 

0 

0 

0 

139 

0 

0 

0 

0 

140 

0 

0 

0 

0 

141 

0 

0 

0 

0 

142 

0 

0 

0 

0 

143 

0 

0 

0 

0 

144 

0 

0 

0 

0 

145 

-0.10198  +  0. 044271  i 

0.11917 -0.098908i 

0.17511  +0.27727i 

-0.17103 -0.30889i 

146 

-0.02217-0.285571 

0.00057447  +  0.516341 

0.023804  +  0.097903i 

-0.01765 -0.1466i 

147 

-0.12407  +  0.1 6402i 

0.69458-0. 36421  i 

0.017724- 0.1 0239i 

-0.047656  +  0.213041 

148 

-2.6184 +  0.14508i 

1.0352- 0.82771  i 

-0.10782  +  0.00369071 

0.52906  +  0.115371 

149 

-0.093182  +  0.15174i 

0.54411  -0.19015i 

-2.995  +  0.173761 

3.5684  +  0.921341 

150 

1.8408  +  1.0087i 

-23.936-  17.67i 

-2.5756- 0.091 033i 

3.9504  +  17.378i 

151 

4.5592  +  0.99307i 

-29.777  -  6.7274i 

-3.0872  +  1.9549i 

17.786  +  3.1931i 

152 

1.2014-6.65441 

-4.6272  +  28.753i 

4.6858  + 2.4081  i 

-10.891  -  14.991i 

153 

0 

0 

0 

0 

154 

0 

0 

0 

0 

155 

-1 

0 

0 

0 

156 

1 

0 

-1 

0 

157 

0 

0 

0 

0 

158 

0 

0 

0 

0 

159 

0 

-1 

0 

0 

160 

0 

1 

0 

-1 
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100 
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B  and  G 

A  and  deIG 

B  and  G 
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APPENDIX  B  -  LEAST  SQUARES  SOLUTION  TECHNIQUE 


The  least  squares  technique  allows  an  overdetermined  system  (“m”  equations  in 
“n”  unknowns,  with  m>n)  to  be  evaluated  in  an  optimized  manner.  This  optimized 
manner  results  in  values  of  x  that  minimize  the  quantity  |Ax  -  b|  for  the  originating 
equation  of  Ax  =  b. 

The  vector  x  that  minimizes  |Ax  -  b|  is  the  solution  to  the  normal  equations: 

A^  Ax  =  A^  b 

T  IT  rvi 

This  vector  x  =  (A*ArA*b  is  the  least  squares  solution  to  Ax  =b  ^  k 

This  appendix  seeks  to  provide  the  reader  with  a  little  more  understanding  of  the 
mechanism  of  this  technique.  A  relatively  simple  example  of  a  sinusoidal  wave 
corrupted  with  a  higher  frequency  random  noise  is  generated.  The  least  squares  solution 
technique  is  then  used  to  determine  a  best  fit  “pure”  sine  wave  to  the  input  data  provided. 
The  sine  wave  so  obtained  minimizes  the  error  per  the  aforementioned  solution 
definitions. 

The  least  squares  solution  technique  can  best  be  envisioned  as  a  means  to  provide 
a  “best”  solution  to  a  long  list  of  equations.  Each  equation  has  the  same  number  and  type 
of  unknown  constants.  For  the  example  chosen,  a  sine  wave  can  be  written  in  terms  of 
both  sine  and  cosine  components.  That  is,  this  signal  can  be  written  as  the  sum  of  these 
two  entities  via: 

D  sin(cot)  +  E  cos(cot)  (1) 

where: 

CO  =  radian  frequency 
t  =  time 
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D  and  E  are  some  defined  eonstants 


In  addition,  the  signal  may  have  a  DC  eomponent,  so  Equation  (1)  may  be  augmented  to: 

D  sin(cot)  +  E  oos(cot)  +  E  (2) 

where: 

E  =  some  eonstant  DC  level 

Equation  (2)  ean  uniquely  speeify  a  given  radian  frequeney  sinusoidally  varying 
signal  by  defining  the  values  for  the  eonstants  D,  E  and  E.  Phase  delays  ean  be 
introdueed  by  the  proper  weighting  of  the  D  and  E  values.  That  is,  if  E=0,  then  the  signal 
so  defined  has  pure  sine  charaeteristies.  If  D  =0,  then  the  signal  so  defined  has  pure 
eosine  eharacteristies.  Non-zero  values  for  both  D  and  E  yield  sinusoidally  varying 
signals  with  some  defined  phase  shift. 

Under  the  assumptions  that  the  radian  frequency,  co,  is  a  known  value  and  that 
time,  t,  can  be  incrementally  adjusted  (ti,  t2,  t3,  U,  etc),  a  set  of  equations  defining  the 
signal  may  be  created. 

D  sin(coti)  +  E  cos(coti)  +  E  (3) 

D  sin(cot2)  +  E  cos(cot2)  +  E 
D  sin(cot3)  +  E  cos(cot3)  +  E 
D  sin(cot4)  +  E  cos(cot4)  +  E 


The  equations  defined  above  consist  of  3  unknowns  (D,  E  and  E)  that  will 
eventually  need  to  be  solved  for.  To  assist  in  this  solution,  the  equations  specified  in  (3) 
above  are  equated  to  a  given  a  set  of  measured  data  values,  m.  This  leads  to: 
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(4) 


D  sin(coti)  +  E  cos(coti)  +  F  =  mi 
D  sin(cot2)  +  E  cos(cot2)  +  F  =  m2 
D  sin(cot3)  +  E  cos(cot3)  +  F  =  m3 
D  sin(cot4)  +  E  cos(cot4)  +  F  =  114 


The  expressions  above  form  the  overdetermined  set  of  equations  that  will  be 
solved  via  least  squares  solution  teehnique.  The  equations  in  (4)  are  rewritten  in  the  form 
of  matriees  yielding: 


sin(coti) 

cos(coti) 

1  *[DEF]  ‘ 

=  mi 

(5) 

sin(cot2) 

cos(cot2) 

1 

=  m2 

sin(cot3) 

cos(cot3) 

1 

=  m3 

sin(cot4) 

cos(cot4) 

1 

=  m4 

where  D  multiplies  all  values  in  column  1  ,  E  multiplies  all  values  in  column  2  and  F 
multiplies  all  values  in  column  3  of  the  first  matrix. 
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The  matrices  in  Equation  (5)  are  now  in  the  form  of  the  defining  equation  Ax  =b  with: 


A  = 

sin(coti) 

cos(coti) 

sin(cot2) 

cos(cot2) 

sin(cot3) 

cos(cot3) 

sin(cot4) 

cos(cot4) 

b  =  mi 

m2 
m3 
rru 


and 

x=  [DEE]’ 

The  solution  may  now  proceed  via  the  least  squares  solution  technique  detailed  at 
the  beginning  of  this  appendix. 

Eor  the  example  chosen,  a  1  Hz  sinusoidally  varying  signal  is  sampled  at  every 
1/256  seconds.  The  data  corresponding  to  the  first  10  rows  (of  the  1024  total)  of  the  A 
matrix  above  was  computed  to  be: 
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A 


1 

0.02454 

0.99970 

1 

2 

0.04907 

0.99880 

1 

3 

0.07357 

0.99729 

1 

4 

0.09802 

0.99518 

1 

5 

0.12241 

0.99248 

1 

6 

0.14673 

0.98918 

1 

7 

0.17096 

0.98528 

1 

8 

0.19509 

0.98079 

1 

9 

0.21910 

0.97570 

1 

10 

0.24298 

0.97003 

1 

The  data  corresponding  to  the  first  10  rows  of  the  b  matrix  were  determined  to  be; 

b 

1  0.68685 

2  0.48398 

3  0.64596 

4  0.64235 

5  0.81405 

6  0.80746 

7  0.74193 

8  0.63214 

9  0.93579 

10  0.84604 


In  like  fashion  as  matrix  A,  vector  b  is  composed  of  1024  total  rows.  The  data 
presented  in  vector  b  consists  of  the  sinusoidally  varying  signal  corrupted  by  some  higher 
frequency  noise.  This  corruption  by  noise  explains  why  the  values  in  b  do  not  vary  in  a 
purely  sinusoidal  manner. 
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The  solution,  x,  for  this  set-up  was  determined  to  be: 


X 


D  1.50520 

E  0.00903 

F  0.50071 


Reeall  that  D  represents  the  sine  eomponent,  E  represents  the  eosine  eomponent 
and  F  represents  the  DC  offset  of  the  measured  signal. 

The  ehart  that  follows  shows  the  original  signal  (blue  crosses)  along  with  the  least 
squares  solution  for  a  best  fit  sinusoidally  varying  signal  (red  line). 


Figure  B1  -  Feast  Squares  Solution  for  Sinewave  with  Noise  and  DC  Offset 


It  is  significant  here  to  address  the  requirement  to  have  an  over  determined  system 
to  utilize  the  method  of  least  squares.  The  sinusoidal  signal  discussed  in  this  section  was 
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represented  by  a  system  of  1024  equations  in  3  unknowns.  This  led,  quite  obviously,  to 
an  over  determined  system  (sinee  m=1024  was  mueh  greater  than  n  =3). 

The  thesis  method  required  the  ineorporation  of  theoretieal  and  physieal 
parameters  to  aehieve  an  over  determined  status.  The  example  afforded  by  the 
AB_Matrix  eonstruetion  of  Appendix  A  ean  best  illustrate  this. 

Appendix  A  deseribed  an  AB_Matrix  with  the  following  eonstituents: 

numv  =  5  ;  this  represents  the  number  of  patehes  in  the  “north-south”  direetion 

sph  sym  =  10  ;  this  represents  the  number  of  “sliees”  that  eomprise  the  sphere 
numh  =  3  ;  this  represents  the  number  of  elements  in  the  near  field  linear  array 

The  eorresponding  system  required  that  a  solution  be  determined  for  50  individual 
patehes  (numv  x  sph  sym).  Eaeh  pateh  required  a  pressure  and  veloeity  solution;  this, 
therefore,  lead  to  the  need  to  solve  for  a  total  of  100  unknowns  (50  pressures  and  50 
veloeities). 


1)  The  use  of  the  first  Helmholtz  relation,  0  =-Ap+Bv,  provides  a  total  of  50 
equations  that  define  the  physieal  system.  The  50  equations  result  from  the 
faet  that  the  AB_Matrix  is  eireulant  and  is  based  upon  the  number  of  vertieal 
patehes  (numv)  and  the  symmetry  order  (sph_sym)  of  the  enelosing  sphere. 
These  50  equations  are,  however,  insuffieient  for  solution  via  least  squares 
sinee  100  unknowns  must  be  solved  for.  This  system  ineorporating  just  this 
Helmholtz  relation  is,  therefore,  under  determined. 

2)  The  use  of  near  field  data  was  essential  to  the  overall  effort  sinee  physieal 
eharaeteristies  of  the  system  needed  to  be  quantified.  The  ineorporation  of 
these  physieal  eharaeteristies  was  embodied  in  the  seeond  Helmholtz  relation, 
pressure  at  a  field  point  =  delG  p  +  G  v.  The  measured  near  field  element 
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data  with  the  above  relationship  formed  another  set  of  defining  equations  for 
the  system. 

The  equations  formed  from  this  relationship  eorresponded  to  the  number  of 
elements  in  the  near  field  array  (numh  =3)  and  the  symmetry  order 
(sph_sym=10).  Thus,  for  the  example,  30  equations  were  assoeiated  with 
this  Helmholtz  relation. 

These  30  equations  are  added  to  the  initial  50  equations  from  the  previous  step 
yielding  80  total  equations  defining  the  system.  These  80  equations,  however, 
still  ean  not  solve  for  the  100  unknowns;  the  system  is  still  under  determined. 

In  this  example,  a  near  field  array  of  3  elements  was  seleeted.  It  would  have 
been  possible  to  define  a  near  field  array  with  more  elements  that  eould  have 
met  the  over  determined  requirement  of  the  least  squares  methodology  (when 
used  with  the  previous  Helmholtz  relations)  .  For  example,  a  10  element 
array  would  have  yielded  100  equations.  These  100  plus  the  previous  50 
would  have  yielded  a  suffieient  number  for  a  least  squares  solution. 

For  this  relatively  simple  example,  a  10  element  array  would  not  have 
presented  any  real  diffieulties  to  implement.  For  other  systems,  the  number  of 
near  field  elements  to  meet  the  over  determined  requirement  eould  approaeh 
staggering  proportions. 

Rather  than  aehieving  the  additional  system  definitions  via  more  near  field 
elements,  a  physies  based  requirement  for  pressures  and  veloeities  on  the 
patehes  was  examined. 

3)  The  pateh  adjaeeney  eonditions  of  pressure  and  veloeity,  pi  =  pi+i  and  Vi  = 
Vi+i,  provided  the  physieal  relationships  to  aehieve  both  additional  defining 
equations  for  the  system  and  to  assist  in  foreing  a  realizable  solution.  Eaeh 
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adjacency  condition  provided  a  total  of  (numv-1)  *  sph  sym  equations.  For 
the  example,  this  yielded  an  additional  40  equations  for  both  pressure  and 
velocity. 


Thus,  the  three  components  above  provided  50,  30,  and  80  (40  and  40)  defining 
equations,  respectively,  for  the  system.  Thus,  an  over  determined  system  was  achieved 
(since  160  equations  now  existed  to  solve  for  100  unknowns). 
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For  completeness  and  further  investigation  ,  the  MATLAB  ending  used  to  define  and 
solve  for  the  sinusoidal  signal  (eorrupted  by  noise)  in  this  appendix  is  provided. 


%funetion  linear_least_squares_fit_to_data_set_SINEWAVE 

%elear  all 
elose  all 

%  Create  Data  Set 
max_points  =  1024; 

time_faetor  =  1/256;  %  sampling  period 

faetor  =3  %  noise  faetor  -  high  numbers  represent  low  noise  -  OPERATOR  INPUT 
radian_phase  =  0;  %  phase  shift  of  eomputed  data  set  -  OPERATOR  INPUT 
offset  =  0.5;  %DC  offset  of  eomputed  data  set  -  OPERATOR  INPUT 
waveform  magnitude  =  1.5;  %magnitude  of  sinewave  -  OPERATOR  INPUT 


for  integer  =I:max_pomts; 
x_value(integer)  =  integer; 
random_number=rand; 

magnitude(mteger)  =  waveform_magnitude.*sm((mteger/256)*2*pi*l+radian_phase)  + 
((random_number-0.5)/faetor)  +  offset;  %+  2*sm((4*integer/256)*pi); 

end 

maxmagnitude  =sqrt(max(magnitude .  ^2)) 
magnitude2  =  magnitude; 

%  Plot  Original  Computed  Data  Set  ineluding  any  noise  eontributions 


142 


figure(l) 

plot  (x_value*time_factor,  magnitude) 


%  Plot  Normalized  Data  Set 
figure(2) 

plot  (x_value*time_faotor,  magnitude!) 
omega  =  2*pi 


time  =[time_faetor:time_faetor:max_pomts*time_factor]  ; 
data  =  magnitude!; 

%  Plot  Normalized  Data  Set  (to  allow  for  direet  eomparison  to  proeessed 

%  data  later) 

figure(3) 

plot(  time, data, 'b+') 
xlabel('time  (t)') 
ylabel('data  value') 

time_size  =  size(time,2) 

%  Set-up  Matriees  and  Solve  for  Least  Square  Coeffieients 
ones_matrix  =  ones([  1  timesize]); 

A  =  [sin(omega*time)  ;  eos(omega*time);  ones  matrix]'; 


B  =  data'; 
A_transpose  =  A'; 
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Result_l  =  A_transpose  *  A; 
Result_2  =  A_transpose  *  B; 
Result_final  =  Result_l  \  Result_2; 


C_value  =  Result_final(l) 

Dvalue  =  Result_final(2) 

Evalue  =  Result_final(3) 

textl  =  ['Least  Squares  Fit  Gives:'] 

text2  =  [  num2str(C_value)  '  sin(wt)  +  '  num2str(D_value)  '  eos(wt)  +  '  num2str(E_value)] 

position_x  =  (max(time)  +  min(time))/8 
position_y  =  (max(data))  %  +  min(data))/2 

%  Plot  data  using  least  squares  derived  eoeffieients 
hold  on 

plot(time,C_value*sin(omega*time)  +  D_value*oos(omega*time)  +  E  value,  'r-') 
text(position_x,  position  y,  text2,  'FontSize',  14) 
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APPENDIX  C  -  EXPLANATION  OF  MATLAB  CODING 


This  appendix  supplements  Chapter  IV  -  Steps  of  the  Methodology.  Detailed 
explanations  of  these  steps  as  implemented  via  the  three  pieees  of  MATLAB  eoding  are 
provided  herein. 

drake_stepl_*,m 

This  step  aeeomplishes  the  following: 

1)  It  ereates  a  simulated  horizontal  transmitting  line  array  of  some  user- 
defined  length.  The  elements  in  the  transmitting  array  are  equally  spaeed  and  are 
spaeed  at  1/6  lambda  maximum.  Lambda  represents  the  wavelength  of  the 
frequeney  of  operation  of  the  array.  This  transmitting  array  definition  would  be 
bypassed  in  the  event  that  an  aetual  deviee  was  used  for  aeoustie  evaluation.  In 
this  latter  ease,  the  maximum  dimension  of  the  aeoustie  deviee  (its  length  for  a 
line  array)  would  be  entered  as  input. 

2)  A  near  field  vertieal  array  is  also  simulated  in  his  seetion  of  the  eode. 
This  is  the  array  that  is  used  in  the  near  field  of  the  aeoustie  deviee  under  test  to 
aequire  response  data.  In  the  eoding  that  is  provided,  the  user  may  define  an  array 
of  arbitrary  length  made  up  of  some  defined  number  of  equally  spaeed  elements. 
In  aetual  praetiee,  the  data  entries  here  would  eorrespond  to  the  aetual 
eonstruetion  details  of  the  near  field  array  being  used. 

3)  An  estimate  of  the  near  field-far  field  boundary  for  the  defined  aeoustie 
deviee  is  ealeulated. 

4)  The  definition  of  range  for  the  plaeement  of  the  near  field  vertieal  array 
for  the  near  field  measurements  is  made.  A  far  field  range  is  also  defined  whieh 
plaees  a  single  reeeive  hydrophone  well  in  the  deviee  under  test’s  far  field.  In  the 
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existing  coding  the  near  field  range  is  defined  as  some  percentage  of  the  distance 
between  the  enclosing  sphere  radius  (described  later)  and  the  near  field-far  field 
boundary  (near  field  range  =  enclosing  sphere  radius  plus  factorA*(nf-ff  boundary 
-  enclosing  sphere  radius);  the  far  field  range  is  represented  as  100  times  this 
latter  factor. 

5)  Test  points  corresponding  to  the  (x,y,z)  positions  for  each  of  the  near  field 
array  vertical  array  elements  rotated  at  some  defined  angular  increment  about  the 
acoustic  device  under  test  are  created.  These  test  points  form  a  cylinder  about  the 
acoustic  device  under  test.  Test  points  for  the  far  field  simulation  are  also  created. 
In  latter  case,  since  the  far  field  will  be  calculated  from  the  response  of  a  single 
hydrophone,  these  test  points  form  a  circle  about  the  device  under  test,  ft  is 
important  to  note  here  that  the  enclosing  sphere  symmetry  order  and  the  number 
of  angular  increments  in  this  section  are  equal. 

In  addition  to  the  test  points  defined  here,  another  set  of  test  points  are 
defined  to  allow  an  evaluation  of  the  overall  technique  when  dimensional 
uncertainty  exists.  The  user  has  the  ability  to  enter  in  values  of  the  error 
associated  with  the  range  to  the  near  field  array,  error  in  its  vertical  placement  and 
errors  in  its  angular  orientation.  In  the  implementation  provided,  these  test  points 
(with  known  errors)  are  associated  with  simulated  data  values  only. 

6)  This  coding  also  calculates  the  complex  pressures  seen  at  each  of  the  near 
field  vertical  array’s  elements  at  each  defined  test  point  and  at  each  of  the  single 
far  field  hydrophone’s  test  points.  The  near  field  test  point  response  data  will  be 
used  to  simulate  actual  measurement  data;  the  far  field  data  so  computed  will  be 
used  later  to  compare  to  the  output  from  the  methodology.  The  calculations  in 
this  step  are  created  by  the  MATLAB  routine  “prsvel”.  This  is  a  custom  function 
call  that  computes  complex  response  given  acoustic  parameters  and  defined 
source  and  field  points.  The  complex  response  at  a  particular  field  point  is 
computed  as  the  contribution  from  each  defined  source  point.  The  complex 
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responses  so  obtained  are  converted  to  dB  values  and  normalized  for  later  display. 
This  section  would  not  be  executed  if  actual  measurement  data  existed. 

7)  The  CHIEF  enclosing  sphere  is  created  within  this  piece  of  coding.  This 
spherical  surface  encloses  the  acoustic  device  under  test  with  a  known 
geometrical  configuration.  This  sphere  will  be  diced  up  into  a  mosaic  of  small 
patches;  a  single  normal  pressure  and  velocity  will  be  assigned  to  each  of  these 
patches  at  a  later  point.  The  radius  of  this  enclosing  sphere  is  deduced  from  the 
size  of  the  acoustic  device  under  test  and  proximity  to  the  location  where  the  near 
field  array  will  be  placed.  In  like  fashion  as  the  placement  of  the  near  field  array, 
the  coding  was  implemented  to  allow  a  percentage  of  the  separation  distance 
between  the  maximum  radial  dimension  of  the  device  under  test  and  the  location 
of  the  near  field  array  to  be  defined  (this  is  factorB  in  the  coding). 

The  number  of  vertical  patches  and  the  enclosing  sphere’s  symmetry  order 
are  also  defined  in  this  section.  The  number  of  vertical  patches  is  based  upon  the 
assumption  that  lambda/6  spacing  is  a  reasonable  maximum  criterion  to  minimize 
pressure  and  velocity  changes  between  adjacent  patches  (this  is  examined  in  more 
detail  below).  The  symmetry  order  is  based  upon  the  number  of  vertical  patches 
with  the  additional  constraint  that  the  symmetry  should  be  a  highly  composite 
number  (a  number  that  may  be  represented  by  powers  of  2  times  3's  and  5's).  This 
latter  requirement  pertains  to  the  EFT  methodology  of  solving  the  matrix 
mathematics  and  will  not  be  addressed  in  detail  here. 

A  plot  of  a  single  section  of  the  enclosing  sphere  is  made  to  visualize  the 
enclosing  surface. 

The  lambda/6  criterion  of  patch  spacing  was  selected  based  upon  an 
empirical  analysis.  This  analysis  looked  at  the  complex  response  at  a  point  in  the 
far  field.  This  response  was  formed  from  the  averaged  sum  of  some  defined 
number  of  common  magnitude  vectors  that  were  equally  spaced  (in  an  arc)  but  all 
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had  different  orientations.  The  resulting  response  in  the  far  field  was  examined 
for  varying  veetor  spaeings  (ie;  varying  fraetional  lambda  spaeings).  This 
response  was  eompared  to  that  which  would  be  achieved  if  all  vectors  were 
perfectly  aligned  with  the  far  field  point  (ie;  all  vectors  pointed  in  the  same 
direction).  This  analysis  showed  that  lambda/6  spacing  yielded  an  error  of  only 
about  0.4  dB  in  the  far  field  response.  This  error  was  deemed  as  satisfactory  for 
the  intended  goals  of  the  technique  and  was  also  expected  to  help  minimize  the 
number  of  patches  required  for  implementation. 

The  figures  below  support  the  above  explanation.  Three  spaeings 
(lambda/6,  lambda/12  and  lambda/4)  are  shown  in  these  figures.  The 
corresponding  difference  in  far  field  magnitude  from  the  first  case  where  all 
vectors  pointed  in  the  same  direction  are  summarized  in  these  figures.  In  these 
examples,  50  vectors  are  shown  in  each  figure. 
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Figure  Cl  -  Vector  Spacing  Summary  and  Effect  on  Far  Field  Pressure 

8)  Using  CHIEF  tools,  the  coding  solves  for  the  A  and  B  matrices  (see 
CHIEF  -  Introduction  for  details  on  these  matrices).  Unlike  previous  versions  of 
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CHIEF,  this  version  allows  a  cessation  of  computations  after  these  matrices  are 
solved  for. 

9)  The  greens  function  integral  function  call  is  next  implemented  for  all 
defined  near  field  test  points  at  a  single  angular  orientation  along  with  all 
enclosing  sphere  patches  (pssph).  This  function  call  results  in  defined  values  of  G 
and  delG  for  later  computations. 

10)  The  Equality  matrix  is  next  established  for  the  pressure  and  velocity 
adjacency  conditions.  This  matrix  consists  has  dimensions  of  (sph  numv-l  x 
sph_sym*sph_numv).  This  arises  because  there  are  sph_numv-l  equality 
conditions  imposed;  each  condition  equates  the  pressure  or  normal  velocity  at  a 
patch  to  its  immediately  adjacent  patch.  This  matrix  takes  the  general  form  of: 


\ — 1 

\ — 1 

1 

0 

0 

0... 

0 

0 

0 

1 

\ — 1 

1 

0 

0... 

0 

0 

0 

0 

0 

0 

0... 

1 

\ — 1 

1 

The  matrix  above  represents  a  section  (sph_numv-l  x  sph  numv)  of  the  overall 
Equality  matrix.  The  remaining  matrix  components  are  all  zeros.  An  Equality 
matrix  for  both  pressure  and  velocity  adjacency  conditions  will  be  provided  later. 

11)  The  final  operation  in  this  piece  of  coding  establishes  a  single  maximum 
value  for  the  A  matrix  and  a  single  maximum  value  from  all  computed  G/rhoc 
and  delG  values.  These  values  will  be  used  in  the  next  step  to  normalize  the 
elements  of  the  combined  AB_Matrix  and  M  matrices. 

12)  The  significant  data  from  this  step  is  saved  to  file. 
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drake_step2_*,m 


This  step  takes  the  data  obtained  from  the  first  step  as  its  input  and  provides  the  following 
funetions: 

1)  The  weighting  factors  alphaO  and  betaO  are  defined  by  the  user.  The 
alphaO  term  is  associated  with  the  theoretical  Helmholtz  Integral  relations  while  the 
betaO  term  is  associated  with  the  near  field  element  data.  (In  practice,  this  latter 
data  set  would  be  as  acquired  by  the  near  field  array;  in  this  effort  the  near  field 
data  set  was  simulated.)  Since  these  values  can  be  selected  by  the  user,  an 
emphasis  can  be  placed  (when  the  pressure  and  velocity  values  of  the  enclosing 
sphere  patches  are  determined)  on  either  the  theoretical  expectations  or  on  the 
empirical  data  obtained.  These  weighting  factors  are  used  along  with  the  values 
from  step  1 1  in  the  first  piece  of  coding  to  create  the  terms  alpha  and  beta. 

2)  The  augmented  Aa  matrix  is  created  by  concatenating  the  A*  alpha, 
delG*beta,  and  Equality  (pressure  and  velocity)  matrices;  the  augmented  Ba  matrix 
is  created  by  concatenating  the  B*alpha/rhoc;  G*beta/rhoc  and  Equality  (pressure 
and  velocity)  matrices.  These  two  matrices  are  then  combined  in  such  a  way  to 
create  the  AB_Matrix  of  form  [AlBl  A2B2A3B3 . . .] 

3)  The  augmented  Ma  matrix  is  next  created  by  concatenating  zeros 
(expected  outputs  for  the  Helmholtz  Integral  relations),  the  near  field  array  element 
data,  and  zeros  (corresponding  to  the  expected  outputs  for  the  (pi  -  pi+i)  and  (vi  - 
Vi+i)  terms). 


4)  The  coding  then  uses  the  “rotational_symmetry_solve”  function  call  to 
create  the  “pv  solution”  matrix.  This  function  call  is  specific  to  the  CHIEE  for 
MATLAB  implementation  and  allows  a  problem  with  rotational  symmetry  to  be 
solved.  Ax  =  B.  A  is  a  (M  x  N)  matrix,  shorthand  for  MxrN  circulant  matrix  where 
r  is  the  order  of  the  symmetry;  b  is  a  (M  by  1)  matrix  of  rhs  values.  The  matrix  so 
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obtained,  pv  solution,  contains  individually  catered  values  of  pressure  and  velocity 
for  each  patch  on  the  enclosing  sphere  (as  determined  by  the  underlying  least 
squares  solution  technique  contained  within  the  “rotational  symmetry  solve” 
function  call.) 

5)  It  is  not  just  sufficient  for  the  solved  for  pressures  and  velocities  to  yield 
values  (when  greens  function  integrations  are  used)  that  agree  with  those  acquired 
by  the  near  field  array  elements.  Without  taking  into  account  the  adjacency 
conditions  and  the  Helmholtz  Integral  relations,  the  solved  for  pressures  and 
velocities  could  be  physically  unrealizable.  To  assess  the  pv  solution  matrix,  the 
function  call  “rotational  symmetry  multiply”  is  used.  This  is  another  function  that 
is  specific  to  the  MATLAB  for  CHIEF  implementation.  This  function  solves  a 
problem  with  rotational  symmetry,  Ax=b.  A  is  a  (M  x  rN)  matrix  where  r  is  the 
order  of  the  symmetry  and  b  is  a  (rN  by  1)  matrix  of  rhs  values.  The  output  from 
this  operation  is  a  new  matrix  designated  as  the  Ma-new  matrix. 

6)  The  balance  of  the  coding  in  this  section  allows  the  operator  to  look  at  the 
constituents  of  both  the  Ma  and  Ma-new  matrices  and  assess  their  agreement. 
Perhaps  the  most  significant  display  in  this  set  pertains  to  the  full  image  outputs  of 
the  Ma  and  Ma-new  matrices.  These  plots  contain  a  lot  of  information  that  allow  a 
direct  comparison  of  the  two  sets  of  receive  array  element  outputs  and  allow  a 
comparison  in  level  of  the  resulting  Helmholtz  Integral  relations  and  adjacency 
conditions  within  the  Ma-new  matrix. 

7)  The  computations  in  this  section  can  be  repeated  with  new  weighting 
factors  (alphaO  and  betaO)  applied  to  try  and  arrive  at  a  best  agreement  between  the 
two  M  matrices.  These  computations  can  be  reiterated  without  the  need  to  re¬ 
create  the  data  set  from  the  first  step. 
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8)  The  primary  output  of  this  coding  is  the  pv  solution  matrix  which  will  be 
used  to  re-establish  pressure  and  velocities  on  the  enclosing  sphere  for  ultimate 
extrapolation  into  the  far  field. 

drake_step3_*,m 

1)  This  coding  takes  the  pv  solution  matrix  from  the  previous  code  and  splits 
it  into  the  pressure  and  velocity  components  corresponding  to  each  of  the 
enclosing  sphere  patches.  Computations  using  greens  functions  integrations  are 
performed  to  arrive  at  the  G  and  delG  values  using  the  enclosing  sphere  geometry 
and  the  mid-point  element  of  the  near  field  array  field  test  points  and  the  defined 
far  field  point.  These  values  of  G  and  delG  are  used  along  with  the  newly 
defiend  pressures  and  velocities  to  compute  pressure  at  a  field  point, 

p  delG  +  V  G. 

2)  Overplots  of  the  original  near  field  data  and  that  obtained  from  this 
technique  are  displayed.  In  addition,  overplots  of  the  original  simulated  far  field 
data  set  and  that  obtained  via  this  technique  are  also  displayed.  In  practice,  this 
last  overplot  would  not  be  made,  instead  just  the  computed  far  field  plot  would  be 
displayed  (since  the  response  would  not  be  known  apriori). 
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%step  #1  -  near  field  to  far  field  extrapolation  techniqne 

close  all 

V.el=25;V.az=25;  %  defines  reference  view  for  3-D  plots 

%  factor  for  placement  of  near  field  array:  sph  rad  -l-factorA(far_field  -  sph  rad) 
factor A=0. 08 

%  factor  for  placement  of  enclosing  sphere:  height/2  -l-factorB(range_near_field  -  height/2) 
factorB=0.1305 

%define  acoustic  media  parameters  and  k/frequency  values  of  operation  for 

%the  transmitting  array  under  test 

ss=1500;rho=1000;rhoc=ss*rho; 

k=20 

freq=k*ss/(2*pi) 
lambda  =  ss/freq 

%define  parameters  for  vertical  near  field  array 

rec  array  spacing  =  [-3*lambda:lambda/2:3*lambda];  %  vertical  spacing  in  meters  of 
%  receive  array  elements 

max  length  rec  array  =  max(rec_array_spacing)  -  min(rec_array_spacing) 
numrecarrayelements  =  size(rec_array_spacing,2) 

%define  parameters  for  horizontal  transmitting  array  (device  under  test) 
height=6*lambda  %  this  is  the  length  of  line  array  in  meters 

num  array  elements  =  (round(height/((l/6)*lambda)))-l-l;  %no.  of  elements  comprising  array 

%  calculate  far  field  criteria  for  defined  line  array  for  info  purposes 
farfieldcomputedcriteria  =  (height^2)/(lambda); 
maxcriteria  =  farfieldcomputedcriteria 

%define  parameters  for  enclosing  sphere  based  upon  wavelength  of 
%transmitting  array  and  its  length 

sph  rad  =  (-(l-factorB)*height/2  -factorA*factorB*max_criteria)/... 
(factorB-(factorA*factorB)- 1) 

%base  number  of  latitude  lines  on  lambda/6  criterion 
sph_numv=ceil(pi*  6  *  sph  rad/ (lambda)) 

%this  gives  a  highly  composite  number  that  speeds  up  FFT  process  later 
sph_sym=2*nearest_small_factor(sph_numv) 

%define  location  of  vertical  near  field  array 

range  near  field  =  -(-factorA*max_criteria  -  sph_rad*(l-factorA)) 

%define  range  where  far  field  is  guaranteed 

range  far  field  =100*  round((max_criteria)/lambda)*lambda  %simulate  far  field 

%  allow  operator  to  assess  number  of  patches  and  other  parameters  before 
%  proceeding 
operator_choice=0 ; 

operator  choice  =  input('Do  you  want  to  use  these  values?  l=yes  0=no  ') 
if  operator  choice  ==  1  %  proceed  with  computations  for  a  "yes"  condition 
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%define  radian  angles  for  which  near  field  data  will  be  required  (the 
%number  of  angles  must  match  the  enclosing  sphere  symmetry) 
th=[0:l:sph_sym-l]'  *(360/(sph_sym))*  pi/180; 

%define  some  degree  of  positional  uncertainty  to  test  methodology 

range  error  =  0;  %  error  in  meters  in  range  to  rec  array 

rec  array  z  error  =  0;  %  error  in  meters  for  vertical  placement  of  rec  array 

th  error  =  0;  %  error  in  radians  of  rec  array  placement 

%define  all  near  field  point  locations 

saved  =  zeros(size(th,l)*num_rec_array_elements,3); 

for  counter=l  :num_rec_array_elements 

test_pts  =  [range_near_field*sin(th),range_near_field*cos(th),... 

zeros(length(th),  1  )+rec_array_spacing(counter)] ; 
test_pts_error  =  [(range_near_field+range_error)*sin(th+th_error)... 
,(range_near_field+range_error)*cos(th+th_error),... 
zeros(length(th),l)+rec_array_spacing(counter)+rec_array_z_error]; 
All_nf_test_pts(  1  +(counter- 1  )*size(th)  :counter*size(th),  1:3)  =test_pts; 
All_nf_test_pts_w_error(  1  +(counter- 1  )* size(th) : counter* size(th),  1:3) ... 
=test_pts_error;  %these  test  points  include  positional  errors 


end 

%defme  all  far  field  point  locations 

test_pts_far=[range_far_field*sin(th),range_far_field*cos(th),zeros(length(th),l)]; 

%  define  all  source  points  of  the  transmitting  line  array 

source_pts  =gen_all_xyz_combinations(0,linspace(-height/2, height/2, num_array_elements),0); 


% 

%surface  definition 

%create  enclosing  sphere  around  array 

pssph=patch_defme('uv_sphere',[-pi/sph_sym,pi/sph_sym,0,pi],l,l,l,sph_numv,... 
[sph_rad,sph_rad,sph_rad], [0,0,0],  [0,0,0]); 

%  p  1  sph  contains  patch  data  for  entire  enclosing  sphere 
p  1  sph=patches_by_rotation(pssph,sph_sym,'z'); 

%  display  single  slice  of  the  enclosing  sphere 
M=plot3_patches(pssph,V); 

%  define  test  points  for  the  entire  enclosing  sphere  based  upon  patch  centers 
[test_pts_encl_sph,norm_encl_sphere,ar_encl_sphere]  =  patch_interp(plsph,.5,.5); 


% - 

%simulate  pressures  at  the  two  defined  ranges  (near  field  and  far  field) 

%define  normal  vector  for  all  defined  source  points 
normals  =  ones(size(source_pts)); 
normals(:,l)  =  0.*normals(:,l); 
normals(:,3)  =  O.*normals(:,2); 

P_prsvel_near_saved  =  []; 

%  compute  pressures  and  velocities  for  all  defined  near  field  test  points 
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for  counter=l  :num_rec_array_elements 

near_field_poiiits  =  All_nf_test_pts_w_error(  1  +(couiiter- 1  )*size(th) : . . . 
counter*  size(th),  1:3); 

[P_prsvel_near,dummy]=prsvel(source_pts,near_field_poiiits,k,rhoc, normals); 
P_prsvel_near_saved(  1  +(counter- 1  )*num_array_elements:counter* . . . 
num_array_elements,l  :size(th))=  P_prsvel_near; 


end 

%  compute  pressures  and  velocities  for  all  defined  far  field  test  points 
[P_prsvel_far,dummy]=prsvel(source_pts,test_pts_far,k,rlioc, normals); 


P_prsvel_near_summed_saved  =  []; 
for  counter=l  :num_rec_array_elements 

summedvalue  =  sum(P_prsvel_near_saved(  1  +(counter- 1 )  *num_array_elements : . . 

counter*num_array_elements,  1 :  size(th,  1 ))); 
P_prsvel_near_summed_saved(counter,:)  =  summedvalue; 
end 

P_prsvel_far_summed  =  sum(P_prsvel_far); 

AlldBPreportednearall  =  []; 

for  counter=l  :num_rec_array_elements 

entry  =  P_prsvel_near_summed_saved(counter,:); 
All_dBP_reported_near_all(counter,:)  =  dB(entry)  -  dB(max(entry)); 
max_rec_array_values(counter,l)  =  dB(max(entry)); 
end 

dBPreportedfar  =  db(P_prsvel_far_summed)-db(max(P_prsvel_far_summed)); 

maxfarfieldvalue  =  db(max(P_prsvel_far_summed)); 

% - 


%Plot  dB  pressures  (simulated)  for  all  near  field  test  points 

for  counter=l  :num_rec_array_elements 
figure 

plot(th*l  80/pi, All_dBP_reported_near_all(counter,:)); 
axis([0  360  -80  0]); 

text2  =  ['NEAR  FIELD  DATA  (Rec  Ele  #'  num2str(counter) ')  with  k  =  '  num2str(k) 
' ;  range  =  '  num2str(range_near_field) '  meters;  array  length  =  ' ... 
num2str(height) '  meters']; 

text3  =  ['  far  field  criteria  =  '  num2str(max_criteria) '  meters;  frequency  =  ' ... 
num2str(freq) '  FIz;  lambda  =  '  num2str(lambda) '  meters'  ]; 

title(text2,  'FontSize',  12) 
text(0,  -10,  text3,  'FontSize',  8) 
xlabel('angle  (degrees)'); 
ylabel('pressure  in  db'); 
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end 


%  Plot  dB  pressures  (simulated)  for  all  far  field  test  points 
figure 

plof(th*  1 80/pi, dBPreporfedfar); 
axis([0  360  -80  0]); 

fexf2  =  ['FAR  FIELD  DATA  wifh  k  =  '  num2sfr(k) ' ;  range  =  '  num2sfr(range_far_field)... 
'  mefers;  array  lengfh  =  '  num2str(height) '  mefers']; 

title(fexf2,  'FontSize',  12); 
fexf(0,  -10,  texfS,  'FontSize',  8); 
xlabel('angle  (degrees)'); 
ylabel('pressure  in  db'); 

%  set  up  CHIEF  parameters  to  eompute  A  and  B  matriees 
elear  O 

0.symmetry='ROTATE'; 

0.symmetry_order=  {spli_sym,'z'} ; 

O .  quadra  ture_order=2 ; 

0.adaptive_quadrature=adaptive_quadrature_strategy(psspli); 

0.rho=rho; 

0.e=ss; 

0.frequeney=freq; 

0.matriees_only=  1 ; 

[p,v]=solve_ehief(psspli,0); 

global  A_  B_ 

N=size(A_,  1 )  ;M=size(A_,2) ; 

%  set  up  to  eompute  greens  funetion  integrals  for  all  defined  near  field 

%  fesf  points 

All_g_nf=[]; 

All_dg_nf=[]; 

for  eounter=l  :num_ree_array_elements 

test_pts  =  All_nf_test_pts(H-(counter-l)*size(th):eounter*size(th),l  :3); 
[g_nf,dg_nf|=greens_funetion_integral(psspli,0,test_pts(l,:)); 

All_g_nf(l:(sph_sym*sph_numv),counter)  =  g_nf; 
All_dg_nf(l:(sph_sym*sph_numv),eounter)  =  dgnf; 

end 

%  Set-up  Equality  Matrix  Strueture 

Ezeros  =  zeros(sph_numv-l,sph_numv); 

Eequality  =  zeros(spli_numv-l,spli_numv); 

for  eounter=l:(sph_numv-l) 

E_equality(eounter,eounter)=l ; 

Eequality  (eounter,eounter-i- 1  )=- 1 ; 
end 

E  temp  =  [E  equality  E  zeros]; 
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for  counter  =  l:sph_sym-2 
Etemp  =  [Etemp  Ezeros]; 
end 

Efmal  =  Etemp; 

Efmalzeros  =  zeros(size(E_fmal)); 

%  Determine  maximum  values  in  A  matrix  and  in  G  and  delG  matrices 
factorA  =  max(max(abs(A_))) 

max_g_nf_templ=  max(abs(All_g_nf(l  :(sph_sym*sph_numv),l))); 
maxdgnftemp  1  =max(abs(All_dg_nf(  1 :  (sphsym*  sphnumv),  1 ))) ; 

for  counter=2 :  (num  rec  array  elements) 

max_g_nf_temp2  =  max(abs(All_g_nf(l  :(sph_sym*sph_numv), counter))); 

max  g  nf  temp  1  =max(max_g_nf_temp  1 ,  max_g_nf_temp2); 

max_dg_nf_temp2=  max(abs(All_dg_nf(l:(sph_sym*sph_numv), counter))); 
max  dg  nf  templ  =max(max_dg_nf_templ,  max_dg_nf_temp2); 

end 

maxgnf  =  abs(max_g_nf_templ); 
maxdgnf  =  abs(max_dg_nf_templ); 

%  use  the  following  factor  to  normalize  the  G  and  delG  data  sets  (note  that 

%  G  also  gets  divided  by  rhoc 

factorgdgnf  =  max(max_dg_nf,max_g_nf/rhoc); 

%clean  up  variables  that  are  no  longer  needed  to  save  memory 

clear  saved  test_ptts  source_pts  near_field_points  dummy  summed  value  P_prsvel_near_saved... 
P_prsvel_far  g_nf  dg  nf  E  temp  entry  E  equality  P_prsvel_near  E  zeros 

%save  the  data  set  to  use  in  the  next  step 

save  drake  full  data  set 

end 
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%step  #2  -  near  field  to  far  field  extrapolation  techniqne 


close  all 

%  load  data  set  from  step  #1 
load  drake  full  data  set 

%  set  to  0  for  first  run,  set  to  any  other  number  for  subsequent  runs  of  step2  to  save  decision  data 
first_run=l; 

if  first_run==0 

saved_step2_data=[]; 

else 

load  step2_decision_data 
end 


%for  group=l  :size(th)  %  set  group  number  for  data  to  be  displayed  -  group  #1 
%  is  the  first  set  of  hydrohone  data  and  its  zeros  range  for  group  is  1  to  size(th) 
group  =1 

%  Set  alphaO  and  betaO  factors  (these  will  be  iterated  for  final  solution) 

%betaO  places  emphasis  on  the  empirical  near  field  data 
betaO  =  1000 

%  alphaO  places  emphasis  on  the  theoretical  Helmholtz  Integral  data 
alphaO  =  1 

%  Compute  final  sclaing  factors  that  will  be  used 
alpha  =  alphaO/factorA 
beta  =  betaO/factorgdgnf 

%create  the  augmented  A  and  B  matrices 

new  A  single  slice  =  [-A_*alpha;All_dg_nf '*beta;  E  final;  E  final  zeros]; 
new  B  single  slice  =  [B_*alpha/rhoc;All_g_nf '*beta/rhoc;  E_final_zeros;E_final]; 

%  interleave  columns  of  new  A  and  new  B  single  slices  into  new  AB  single  slice 
ABMatrix  =  [newAsingleslice  newBsingleslice]; 
clear  new  A  single  slice  new  B  single  slice 

%  now,  put  indices  in  order  to  keep  circulant: 

%  interleave  from  [A1A2A3...B1B2B3...]  to  [A1B1A2B2A3B3...]  (block 
%  circulant  form),  where  each  Am  amd  Bn  is  numv  by  numv  in  size, 
perm_ix=reshape([l  :M;M+[1  :M]],2*M,1); 

AB_Matrix=AB_Matrix( :  ,perm_ix) ; 


%  Set  up  rhs  of  equation  ,  Ma,  the  vector  of  actual  measurements 

%  Define  two  vectors  of  zeros  corresponding  to  A  and  B  (number  of  rows)  and 
%  the  Equality  Matrix  (number  of  rows) 

ZeroO  =  zeros(sph_numv,l); 

Zerol  =  zeros(sph_numv-l,l); 

Mtemp  =  []; 

%  build  up  M  matrix  using  ZeroO,  Zerol  and  the  measured  near  field 

158 


%  pressures  at  eaeh  reeeive  element 


for  eounter  =l:size(th); 

%  "P_prsvel_near_summed"  represents  the  eomplex  pressure  measured  at  a 
%  partieular  field  point.  In  this  eode,  it  was  eomputed  via  "prsvel") 

M  temp  =  [M  temp;  ZeroO;P_prsvel_near_summed_saved(:, eounter);  Zerol;  Zerol]; 
end 

Mfinal  =  M_temp*beta; 
elear  Mtemp 

group  size  =  (3*sph_numv)  -2  +  size(ree_array_spacing,2); 

%solve  for  pressure  and  veloeity 

pv  solution  =  rotational_symmetry_solve(AB_Matrix,  M  final,  sph  sym); 

%eompute  Ma-new  using  solved  for  pressures  and  veloeities 

M_from_pv_soln=rotational_symmetry_multiply(AB_Matrix,pv_solution,sph_sym); 

%  separate  eonstituents  from  the  two  M  matriees 
M  final  resizel  =  reshape(M_final,  group  size,  sph  sym); 

M_final_resize2  =  M_final_resizel(sph_numv+1  :sph_numv+size(ree_array_spaeing,2),:); 

M_final_resize3  =  reshape(M_final_resize2,(size(ree_array_spaeing,2)*sph_sym),  1); 

M_final_zero0_resize2  =  M_fmal_resizel(l:sph_numv,:); 

M_final_zero0_resize3  =  reshape(M_final_zeroO_resize2,(sph_numv*sph_sym),  1); 

Mfinalzero  1  _resize2  = 

Mfmalresize  1  (sph_numv+size(ree_array_spaeing,2)+l  :3  *sph_numv+size(ree_array_spaeing,2)-2, :); 
M_final_zerol_resize3  =  reshape(M_final_zerol_resize2,((2*sph_numv-2)*sph_sym),  1); 

M_from_pv_soln_resizel  =  reshape(M_from_pv_soln,  group  size,  sph  sym); 

M_from_pv_soln_resize2  = 

M_from_pv_soln_resize  1  (sph_numv+ 1 :  sph_numv+size(ree_array_spaeing,2), : ) ; 

M_from_pv_soln_resize3  =  reshape(M_from_pv_soln_resize2,(size(ree_array_spaeing,2)*sph_sym),  1); 

M_from_pv_soln_zero0_resize2  =  M_from_pv_soln_resizel(l  :sph_numv,:); 

M_from_pv_soln_zero0_resize3  =  reshape(M_from_pv_soln_zeroO_resize2,(sph_numv*sph_sym),l); 

M_from_pv_soln_zerol_resize2  = 

M_from_pv_soln_resizel(sph_numv+size(ree_array_spaeing,2)+l:3*sph_numv+size(ree_array_spaeing,2) 

-2,:); 

M_from_pv_soln_zerol_resize3  =  reshape(M_from_pv_soln_zerol_resize2,((2*sph_numv- 

2)*sph_sym),l); 

maxhydro  =  10*(round((max(dB(abs(M_final_resize3))))/10)) 

ymax  =  maxhydro; 
ymin  =  max  hydro  -  150; 

figure 

plot(l:size(ree_array_spaeing,2)*sph_sym,dB(abs(M_final_resize3)),'k.-', 
l:size(ree_array_spaeing,2)*sph_sym,dB(abs(M_from_pv_soln_resize3)),'r.-' ) 

textlO  =  ['ALL  HYDROPHONE  DATA  (dB);  alphaO  =  '  num2str(alpha0)  '  ;  alpha  =  '  num2str(alpha)  '  ; 
betaO  =  '  num2str(beta0) ' ;  beta  =  '  num2str(beta)]; 
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title(textlO) 

axis([l  size(rec_array_spacing,2)*sph_sym  ymin  ymax]) 
xlabel('Near  Field  Array  Field  Points') 
ylabel('Level  (dB)') 

figure 

plot(th*  1 80/pi, min(dB(abs(M_from_pv_soln_resize2))),'r.-' ) 

textlO  =  ['min(Ma-new  FIYDROPFIONE  DATA)  (dB);  alphaO  =  '  num2str(alpha0)  '  ;  alpha  = 

num2str(alpha) ' ;  betaO  =  '  num2str(beta0) ' ;  beta  =  '  num2str(beta)]; 

title(textlO) 

axis([0  400  ymin  ymax]) 
xlabel('Angle  (degrees)') 
ylabel('Level  (dB)') 

M  hydro  difference  =  dB((M_from_pv_soln_resize3  -  ]Vt_final_resize3)/max(max(]Vt_final_resize3))); 
M_hydro_difference2  =  dB(]Vt_from_pv_soln_resize3)  -  dB(]Vt_final_resize3); 

figure 

plot(l:size(rec_array_spacing,2)*sph_sym,M_hydro_difference,'g.-'); 

textlO  =  ['HYDROPFIONE  DATA  -  dB[(computed-actual)/max(actual)]  ;  alphaO  =  '  num2str(alpha0)  ' 

betaO  = '  num2str(beta0)]; 

title(textlO) 

xlabel('Near  Field  Array  Field  Points') 
ylabel('Level  (dB)') 


figure 

plot(l:size(rec_array_spacing,2)*sph_sym,M_hydro_difference2,'g.-'); 

textlO  =  ['FIYDROPFIONE  DATA  -  dB[(computed-actual)]  ;  alphaO  =  '  num2str(alpha0)  '  ;  betaO  = 

num2str(beta0)]; 

title(textlO) 

xlabel('Near  Field  Array  Field  Points') 
ylabel('Level  (dB)') 


figure 

plot(l:sph_numv*sph_sym,dB(abs(M_from_pv_soln_zero0_resize3)),'b.-') 

textlO  =  ['ALL  Ma-new  ZeroO  Data  (Helmholtz  Integral)  (dB);  alphaO  =  '  num2str(alpha0)  '  ;  betaO  = 

num2str(beta0)]; 

title(textlO) 

axis([l  sph_numv*sph_sym  ymin  ymax]) 
ylabel('Level  (dB)') 

xlabel('Helmholtz  Integral  Conditions  (sphnumv  x  sphsym)') 
figure 

plot(  1 :  (2  *  sph_numv-2)  *  sph_sym,dB(abs(M_from_pv_soln_zero  1  _resize3 )),  'k.  -') 

textlO  =  ['ALL  Ma-new  Zerol  DATA  (Pressure/Velocity  Equivalence)  (dB);  alphaO  =  '  num2str(alpha0) ' 

betaO  = '  num2str(beta0)]; 

title(textlO) 

axis([l  (2*sph_numv-2)*sph_sym  ymin  ymax]) 
ylabel('Level  (dB)') 

xlabel('Pressure/Velocity  Adjacencies  ((2*sphnumv-2)  x  sphsym)') 
figure 

plot(th*  1 80/pi, max(dB(abs(Mfrom_pv_soln_zeroO_resize2))),'b.-') 

160 


textlO  =  ['max(Ma-new  ZeroO  Data)-  (Helmholtz  Integral)  (dB);  alphaO  =  '  num2str(alpha0)  '  ;  betaO  =  ' 

num2str(beta0)]; 

title(textlO) 

axis([0  400  ymin  ymax]) 
xlabel('Angle  (degrees)') 
ylabel('Level  (dB)') 

figure 

plot(th*  1 80/pi, max(dB(abs(M_from_pv_soln_zero  l_resize2))),'k.-') 

textlO  =  ['max(Ma-new  Zerol  DATA)  -  (Pressure/Velocity  Equivalence)  (dB);  alphaO  =  '  num2sb(alpha0) ' 

;  betaO  =  '  num2str(beta0)]; 

title(textlO) 

axis([0  400  ymin  ymax]) 
xlabel('Angle  (degrees)') 
ylabel('Level  (dB)') 

M  hydro  zeroO  diff  =  dB(min(M_from_pv_soln_resize2))  -  dB(max(M_from_pv_soln_zeroO_resize2)); 

M  hydro  zerol  diff  =  dB(min(M_from_pv_soln_resize2))  -  dB(max(M_from_pv_soln_zerol_resize2)); 

figure 

plot(th*  1 80/pi, M_hydro_zeroO_diff,'g.-'); 

textlO  =  ['Ma-new  [  min(HYDROPHONE)  -  max(ZeroO)  ]  (dB)  ;  alphaO  =  '  num2str(alpha0)  '  ;  betaO  =  ' 

num2sb(beta0)]; 

title(textlO) 

xlabel('Angle  (degrees)') 
ylabel('Level  (dB)') 

figure 

plot(th*  1 80/pi,M_hydro_zero  l_diff,'g.-'); 

textlO  =  ['Ma-new  [  min(HYDROPHONE)  -  max(Zerol)  ]  (dB)  ;  alphaO  =  '  num2str(alpha0)  '  ;  betaO  =  ' 

num2sb(beta0)]; 

title(textlO) 

xlabel('Angle  (degrees)') 
ylabel('Level  (dB)') 


figure 

plot(th*  1 80/pi,min(dB(abs(M_from_pv_soln_resize2))),'r.-',... 
th*  1 80/pi,max(dB(abs(M_from_pv_soln_zero0_resize2))),'b. 
th*180/pi,max(dB(abs(M_from_pv_soln_zerol_resize2))),'k.-') 
textlO  =  ['Ma-new  -  min(HYDROPHONE),  max(ZeroO),  max(Zerol);  alphaO  =  '  num2sb(alpha0)  '  ;  alpha 
=  '  num2str(alpha) ' ;  betaO  =  '  num2str(beta0) ' ;  beta  =  '  num2sb(beta)]; 
textll  =  ['  red:  min(hydrophone);  blue:  max(ZeroO);  black:  max(Zerol)']; 

title(textlO) 

axis([0  400  ymin  ymax]) 
xlabel('Angle  (degrees)') 
ylabel('Level  (dB)') 
text(0,  ymin-i-10,  textll,  'FontSize',  8) 

dim  =[ymin,ymax] 
figure 

subplot(l,2,l);imagesc(max(-200,db(M_from_pv_soln_resizel)),cbm) 

colorbar 

title('Ma-new  -  from  pv  solution') 
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subplot(l,2,2);imagesc(max(-200,db(M_final_resizel)),clim) 

colorbar 

title('Ma  -  original') 

%clean  up  variables  that  are  no  longer  needed 

clear  A_  B_  All  dg  nf  All  g  nf  E  final  E  fmal  zeros  M  final  M_from_pv_soln 


%  compute  parameters  to  be  saved  and  examined  later 

savedl  =  min(min(dB(abs(M_from_pv_soln_resize2))))  %minimum  of  the  new  hydrophone  data  minimum 
set 

saved2  =  max(min(dB(abs(M_from_pv_soln_resize2))))  %  maximum  of  the  new  hydrophone  data 
minimum  set 

savedS  =  mean(min(dB(abs(M_from_pv_soln_resize2))))  %average  of  the  new  hydrophone  data  minimum 
set 

saved4  =  min(max(dB(abs(M_from_pv_soln_zeroO_resize2))))  %minimum  of  the  new  zeroO  data 

maximized  set 

saved5  =  max(max(dB(abs(M_from_pv_soln_zeroO_resize2))))  %maximum  of  the  new  zeroO  data 
maximized  set 

savedb  =  mean(max(dB(abs(M_from_pv_soln_zeroO_resize2))))  %average  of  the  new  zeroO  data 

maximized  set 

saved?  =  min(max(dB(abs(M_from_pv_soln_zerol_resize2))))  %minimum  of  the  new  zerol  data 

maximized  set 

savedS  =  max(max(dB(abs(M_from_pv_soln_zerol_resize2))))  %maximum  of  the  new  zerol  data 
maximized  set 

saved9  =  mean(max(dB(abs(M_from_pv_soln_zerol_resize2))))  %average  of  the  new  zerol  data 

maximized  set 

savedlO  =  min(M_hydro_difference2)  %minimum  of  the  (new-original)  hydrophone  data  set 
savedl  1  =  max(M_hydro_difference2)  %maximum  of  the  (new-original)  hydrophone  data  set 
savedl2  =  mean(M_hydro_difference2)  %average  of  the  (new-original)  hydrophone  data  set 

savedl3  =  min(M_hydro_zeroO_diff);  %minimum  of  the  (hydrophone(new)  -  zeroO(new)  data 
savedl4  =  max(M_hydro_zeroO_diff);  %maximum  of  the  (hydrophone (new)  -  zeroO(new)  data 

saved  15  =  min(M_hydro_zerol_diff);  %minimum  of  the  (hydrophone(new)  -  zerol  (new)  data 
savedlb  =  max(M_hydro_zerol_diff);  %maximum  of  the  (hydrophone (new)  -  zerol(new)  data 


%  order  of  data  - 

%  alphaO  (1) 

%betaO  (2) 

%  min  hydro  diff  (3) 

%  max  hydro  diff  (4) 

%  mean  hydro  diff  (5) 

%  min  hydro  minus  zeroO  (6) 
%  max  hydro  minus  zeroO  (7) 
%  min  hydro  minus  zero  1  (8) 
%  max  hydro  minus  zero  1  (9) 
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%  min  new  hydro  set  (10) 

%  max  new  hydro  set  (11) 

%  mean_new_hydro_set  (12) 

%  min  new  zeroO  set  (13) 

%  max_new_zeroO_set  (14) 

%  mean_new_zero0_set(15) 

%  min  new  zerol  set  (16) 

%  max_new_zerol_set  (17) 

%  mean_new_zerol_set  (18) 

step2_data  =  [alphaO  betaO  saved  10  saved  11  saved  12  saved  13  saved  14  ... 
saved  15  saved  16  saved  1  saved2  saved3  saved4  saved5  saved6  ... 
saved?  savedS  saved9]; 

if  first  run  ==0 

saved_step2_data  =  [step2_data] 
else 

saved_step2_data  =  [saved_step2_data;  step2_data]; 
end 

header_info={'alphaO  (l)VbetaO  (2)  Vmin  hydro  diff  (3)',  'max  hydro  diff  (4)Vmean_hydro_diff  (5)',... 
'min  hyd  -  zeroO  (6)','  max  hyd  -  zeroO  (7)Vinin_hyd_-_zerol  (8)'  'max  hyd  -  zerol  (9)'} 

header=(l:9); 

[header;  saved_step2_data(:,  1:9)] 


save  step2_decision_data  saved_step2_data; 
wklwrite('step2_data_wkl_fde',  saved_step2_data) 
save  drake  full  data  set2 
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%step  #3  -  near  field  to  far  field  extrapolation  techniqne 


close  all 

%load  data  set  from  step  #2 
load  drake  full  data  set2 


%  pull  pressure  and  velocity  values  from  previous  step  for  use  in 
%  green's  function  integrations 

pv_solnl=  resliape(pv_solution,2,size(pv_solution,l)/2); 

penclsphere  =  pv_solnl(l,:).'; 
venclsphere  =  pv_solnl(2,:).'/rhoc; 

clear  pv  solnl  pv  solution 


%select  near  field  array  element  in  the  middle 
counter  =  ceil(num_rec_array_elements/2) 

%evaluate  green's  function  from  all  patches  to  all  near  field  points  and 
%calculate  pressures  for  all  test  points  for  the  middle  recieve  array  element 

test_pts  =  All_nf_test_pts(l+(counter-l)*size(th,l):counter*size(th,l),l  :3); 

[g_nf,dg_nf|=greens_fiinction_integral(pssph,0,test_pts); 

fp_temp=calculate_field_pressures(dg_nf,g_nf,p_encl_sphere,v_encl_sphere); 

max_field_point  =  max(fp_temp); 

Field_pressures(counter,l:size(th,l))  =  dB(fp_temp)  -  max_rec_array_values(counter,l); 
clear  fp  temp  g_nf  dg  nf 

%Plot  dB  Pressures  for  near  field  (original  and  that  obtained  via 
%technique) 

figure 

plot(th*l 80/pi, All_dBP_reported_near_all(counter,l:size(th,l)),'b.-',  ... 
th*  1 80/pi, Field_pressures(counter,  1 :  size(th,  1 )),  'r. : '); 

axis([0  360  -80  0]); 

text2  =  ['NEAR  FIELD  DATA  (Rec  Ele  #'  num2str(counter) ')  with  k  =  ' ... 

num2str(k) ' ;  range  =  '  num2sfr(range_near_field) '  meters;  array  length  =  ' ... 
num2str(height) '  meters']; 

text3  =  ['  far  field  criteria  =  '  num2sfr(max_criteria) '  meters;  frequency  =  ' ... 

num2str(freq) '  FIz;  lambda  =  '  num2sfr(lambda) '  meters'  ]; 
text4  =  ['  alpha  =  '  num2sfr (alpha) ' ;  beta  =  '  num2str(beta) ' ;  alphaO  =  ' ... 
num2str(alpha0) ' ;  betaO  =  '  num2str(beta0)  ] 

text6  =  ['  number  of  equally  spaced  receive  array  elements:  '  num2str(num_rec_array_elements)... 

' ;  its  length  is  '  num2sfr(max_length_rec_array) '  meters']; 
texts  =  ['  red  (computed  via  method) ;  blue  (measured/simulated)']; 
text?  =  ['  enclosing  sphere  radius  =  '  num2sfr(sph_rad) '  meters;  sph  numv  =  ' ... 
num2str(sph_numv) ' ;  sph  sym  =  '  num2sfr(sph_sym)]; 
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title(text2,  'FontSize',  12) 
text(0,  -10,  texts,  'FontSize',  8) 
text(0,  -70,  text4,  'FontSize',  8) 
text(0,  -60,  texts,  'FontSize',  8) 
text(0,  -65,  text6,  'FontSize',  8) 
text(0,  -75,  text7,  'FontSize',  8) 

xlabel('angle  (degrees)'); 
ylabel('pressure  in  db'); 


%Plot  dB  Pressures  differenees  for  near  field  (original  and  that  obtained  via 
%teehnique) 

figure 

plot(th*180/pi,All_dBP_reported_near_all(eounter,l:size(th,l))-  Field_pressures(eounter,l:size(th,l)),'g.-'); 

text2  =  ['NF  DATA  Differenee  (Ree  Ele  #'  num2str(eounter) ')  with  k  =  ' 
num2str(k) ' ;  range  =  '  num2str(range_near_field) '  meters']; 

title(text2,  'FontSize',  12) 

xlabel('angle  (degrees)'); 
ylabel('db'); 


%evaluate  green's  funetion  and  ealeulate  pressure  for  all  far  field  poinfs 

[g_ff,dg_ff|=greens_funefion_integral(pssph,0,fesf_pfs_far); 

P_ff_gfl=ealeulafe_field_pressures(dg_ff,g_ff,p_enel_sphere,v_enel_sphere); 

%normalize  far  field  pressures  eompufed  via  Greens  Funetion  and  eonvert  to  dB 
%P_est_reported_ffl  =  db(P_ff_gfl)-db(max(P_ff_gfl)); 

Pestreportedffl  =  db(P_ff_gfl)-max_far_field_value; 

elear  g_ff  dgff 


%Plot  dB  Pressures  for  far  field  (original  and  that  obtained  via 

%teehnique) 

figure 

plot(th*  1 80/pi,dBP_reported_far,'b.-',th*  1 80/pi,P_est_reported_ffl ,  'r. :' ); 
axis([0  360  -80  0]); 

text2  =  ['FAR  FIELD  DATA  (on  axis)  with  k  =  '  num2str(k) ' ;  range  =  ' ... 

num2str(range_far_field) '  meters;  array  length  =  '  num2str(height) '  meters']; 
texts  =  ['  far  field  eriteria  =  '  num2str(max_eriteria) '  meters;  frequeney  =  ' ... 

num2str(freq) '  Flz;  lambda  =  '  num2str(lambda) '  meters'  ]; 
text4  =  ['  alpha  =  '  num2str(alpha) ' ;  beta  =  '  num2str(beta) ' ;  alphaO  =  ' ... 

num2str(alpha0) ' ;  betaO  =  '  num2str(beta0)  ] 
text5  =  ['  red  (eomputed  via  method) ;  blue  (measured/simulated)']; 

text6  =  ['  number  of  equally  spaeed  reeeive  array  elements:  '  num2str(num_ree_array_elements)... 

' ;  its  length  is  '  num2str(max_length_ree_array) '  meters']; 
text7  =  ['  enelosing  sphere  radius  =  '  num2str(sph_rad) '  meters;  sph  numv  =  ' ... 
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num2str(sph_numv) ' ;  sph  sym  =  '  num2str(sph_sym)]; 


title(text2,  'FontSize',  12) 
text(0,  -10,  texts,  'FontSize',  8) 
text(0,  -70,  text4,  'FontSize',  8) 
text(0,  -60,  texts,  'FontSize',  8) 
text(0,  -65,  text6,  'FontSize',  8) 
text(0,  -75,  text7,  'FontSize',  8) 

xlabel('angle  (degrees)'); 
ylabel('pressure  in  db'); 


%Plot  dB  Pressures  differenees  for  far  field  (original  and  that  obtained  via 

%teehnique) 

figure 

plot(th*180/pi,dBP_reported_far-P_est_reported_ffl,  'g.-' ); 

text2  =  ['FF  DATA  Differenee(on  axis)  with  k  =  '  num2str(k) ' ;  range  =  ' 
num2str(range_far_field) '  meters']; 

title(text2,  'FontSize',  12) 

xlabel('angle  (degrees)'); 
ylabel('db'); 


%  Re-plot  All  Flydrophone  Data  that  eorresponds  to  the 
%  data  set  presented. 

figure 

plot(l:size(ree_array_spaeing,2)*sph_sym,abs(M_final_resize3), ... 
l:size(ree_array_spaeing,2)*sph_sym,abs(M_from_pv_soln_resize3) ) 

textlO  =  ['ALL  HYDROPHONE  DATA  ;  alphaO  =  '  num2str(alpha0) ' ;  alpha 
num2str(alpha) ' ;  betaO  =  '  num2str(beta0) ' ;  beta  =  '  num2str(beta)]; 
title(textlO) 


elim  =[ymin,ymax] 
figure 

subplot(l,2,l);imagese(max(-200,db(M_from_pv_soln_resizel)),ehm) 

eolorbar 

title('Ma-new  -  from  pv  solution') 

subplot(l,2,2);imagese(max(-200,db(M_final_resizel)),ehm) 

eolorbar 

title('Ma  -  original') 

range  error  %  error  in  meters 
ree  array  z  error  %  error  in  meters 
th  error  %  error  in  radians 

save  drake  full  data  set  final 
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